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te. INTRODUCTION 


Compressor stall and off-design behavior limit the per- 
formance attainable in the design of aircraft gas turbine 
engines. The occurrence of compressor stall in flight can 
lead to a loss of engine power and, in some circumstances, 
can lead to the loss of the rare The need to develop 
higher thrust-to-weight ratio engines for  fiahter aircraft, 
and distortion tolerant engines for vertical landinq 
applications, requires that analytical tools be developed to 
describe the phenomenon of stall, and particularly, its 
onset. Ideally the compressor flow field should be fully 
predictable and quantifiable in the design process at all 
possible operating conditions. 

The aerodynamic design of compressor cascades, which is 
equivalent to designing the blade elements on axisymmetric 
stream surfaces through the machine, is usually based on 
simplified but realistic modelling, implemented using 
Computational Fluid Dynamic (CFD) codes. Desian point is 
analyzed using Euler or potential methods, addina a boundary 
layer prediction, and usina empirical separation criteria to 
avoid stall. The prediction of the off-desian behavior and 
stall onset however, requires much more sophisticated CFD 
methods. Such predictions are clearly essential in the 


selection of optimum designs. Viscous codes based on the 


full Navier-Stokes (N-S) equations are becoming increasingly 
available, but are generally too əxpensive for routine 
application in desiqn refinement. In order to overcome the 
ә of the N-S codes, but to obtain a reasonably 
accurate description of the flov field containing small 
regions of separated flow, or simply thick boundary layers, 
viscous-inviscid strong interaction codes are currently 
being developed which incorporate either integral or differ- 
ential boundary layer descriptions. Unfortunately, all com- 
putational descriptions of viscous compressor cascade flow 
fields must rely on empirical models for transition and 
turbulence. Clearly, the CFD code will only predict the 
cascade flow well if the empirical models are qood ones. 
Thus, what is critical to obtaining accurate predictions 
of the off-design behavior of compressor cascades are care- 
fully controlled experiments in which the viscous code 
descriptions are, in effect, calibrated. It was found that 
there were no data available in the open literature from 
realistic compressor blade geometries which were in suffi- 
cient detail to be used for code validation or calibration. 
Therefore, the present study, in vhich the emphasis vas 
placed first on obtaining the needed data, was initiated. 
What was achieved in the study was 
a. a detailed mapping of the flow field and, in precise 


detail, the viscous layer development ina modern 


controlled-diffusion compressor cascade at design and 
two higher incidence angle, off-design conditions, 
a comparison of viscous layer growth with predictions 
obtained using one generally available, non- 
interactive code, a fully interactive code, and a 
Navier-Stokes code; 

the identification of limitations and critical needs 


common to all CFD approaches. 


II. NATURE OFE THE PROBEEM 


A. CONTROLLED-DIFFUSION (CD) BLADING DESIGN 

CD bladinqs have been developed in the last decade to 
replace the standard series airfoils such as NACA-65 series, 
Double Circular Are (DCA), or Multiple Circular Arc (CEE 
The CD design vas enabled by the evolution of CPD 
techniques. Such methods, initially based on inviscid flow 
analysis, could be used to arrive at blade surface shapes 
having prescribed pressure distributions. By controlling 
diffusion on the suction surface of the airfoil, boundary 
layer separation could be avoided at design point, and the 
flow could be designed to be shock-free in che transonic 
range. 

CD blades are quite different from NACA-65 series or DCA 
airfoils in that the shapes of the upper and lower surfaces 
are "arbitrary," resulting from the use of high order 
polynomials. Also, CD airfoils can be considerably more ro- 
bust toward the leading and trailing edges, especially com- 
pared to the DCA blading. As a result, the trailing edae 
radius may also be larger in the CD airfoil. 

Several organizations have developed methods to desiqn 
CD airfoils but have reported the results only recently. 
Dunker et al (Ref. 1) designed and tested CD blade shapes in 


a redesigned compressor stator and in a plane cascade wind 


tunnel at the DFVLR in West Germany. In the U.S., Hobbs and 
Weingold (Ref. 2) developed a CD blade design system for 
multistage compressor applications at Pratt & Whitney 
2 craft. The design vas based on boundary layer calcula- 
tions incorporating laminar-turbulent transition by Dunham 
(Ref. 3) and a laminar separation and reattachment criterion 
by Roberts (Ref. 4). At NASA's Lewis Research Center, 
Sanger (Ref. 5) used an optimization technique to design CD 
cascades and applied the method in the redesign of the first 
Stage stator of a two-stage fan. Sanger used a series of 
computational methods coupled by a numerical optimization 
procedure. The blade shape, from the mean line shape and 
thickness distribution, was defined using high order 
polynomials. Potential flow was calculated by a code devel- 
oped by Katsanis (Ref. 6) and boundary-layer calculations 
were carried out using BLAYER by McNally (Ref. 7). The CD 
test blades and cascade geometry used in the present study 
were scaled from the midspan section of Sanger's redesigned 
stator. 

All published design methods have used an inviscid flow 
calculation to predict pressure distribution followed by a 
boundary layer calculation to predict viscous behavior. The 
surface contour is generally modified to account for bound- 
ary layer displacement thickness anda criterion of some 
kind is adopted to allow designs to approach but not to 


exceed the suction side adverse pressure gradient which 


vould result in separation. For as long as separation is 
avoided, and the boundary layers are thin compared co the 
blade spacing, such modelling is appropriate, and gives good 
results. At well off-design conditions, corresponding to 
reduced through-flow, thickening boundary layers and the on- 
set of separation make such modelling no longer accurate. 
To date reliable methods for the prediction of the off- 
design and, particularly, near-stall behavior of CD blading 


do not exist. 


B. DESCRIPTION OF THE OF F-DEStGC T STT TETE 

A schematic to illustrate features of the flow field 
over a CD blade at an incidence angle within its working 
range is shown in Figure l. The significant features are 
labelled from 1 to 8 as follows: 

1. The stagnation point, normally within the leading edge 
curvature. 

2. Rapid acceleration of the flow around the leading edge 
radius. The boundary layer here is laminar. 

3. Laminar separation, creating a separation bubble. 
Transition occurs in the free shear layer with re- 
attachment usually as a turbulent boundary layer. An 
enlargement of this region is shown as described by 
Walker (Ref. 8). 

4. Deceleration of the flow in an adverse pressure gradi- 


ent with a turbulent boundary layer. 
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5. Separation of the turbulent boundary aee 

SA. When the turbulent boundary-layer remains attached 
over the suction surface, SIM occurs at the 
curvature of the blunt trailing edge, creating back= 
flows in the near wake. 

5B. Mixing in the far wake. 

6. Acceleration of the flow on the leading edge radius of 
the pressure surface followed by a fairly uniform 
velocity distriDurfiom: The boundary layer here is 
normally laminar. 

7. ‘Boundary layer transition. 

Sa TurDu len boundary layer under a fairly uniform 
velocity distriDutTon” 

C. CALCULATING THE OFF=|DESIGN FLOW BEHAVIOR 

During the past twenty-five years, computational tech- 
niques for computing flow fields have become increasingly 
more sophisticated. Inviscid codes are well developed and 
used extensively in the design process. Viscous codes based 
on the boundary layer equations are used when the boundary 
layer is attached and relatively thin and the outer, invis- 
cid flow has been solved. Such an approach is also appro- 
priate for conditions not far from the desiqn point if addi- 
tional modelling is used for the wake flow. The viscous 
code BLAYER was developed for the analysis of cascade flow 
by McNally (Ref. 7). The BLAYER transition is predicted 


using the Schlichting-Ulrich-Granville method (Ref. 9). 


Laminar separation creatinq a separation bubble is assumed 
to occur when wall shear stress passes from positive to 
negative. Turbulent separation is assumed to occur when the 
magnitude of the incompressible shape factor reaches a 
specified value. 

Calculations based on the solution of the Navier-Stokes 
equations can, in principle, be used to calculate both on- 
and off-design conditions and to model, numerically, all the 
features of the flow field shown in Figure l. The litera- 
ture includes papers by Beam and Warming (Ref. 10), Steger 
Meet. 11), Rubin and Koshla (Ref. 12) to solve the external 
flow around airfoils and Shamroth, McDonald, and Briley 
(Ref. 13), and Shafer et al (Ref. 14) for internal flow and 
turbomachinery applications. 

To overcome the extended computer time and memory re- 
quirements of Navier-Stokes codes calculation techniques 
based on  viscous-inviscid interactive procedures have been 
developed for internal flows in recent years. Johnson and 
Sockol (Ref. 15) developed an interactive code for a cascade 
el airfoils. Cebeci, Hess, and Lee (Ref. 16) adopted an 
external aerodynamic interactive code for cascade applica- 
tions, and Krainer (Ref. 17) extended the code to calculate 
the flow in a NACA-65 series airfoil cascade. Viscous- 
inviscid interaction techniques appear to work when regions 
of separation are small, and therefore, should be useful in 


the prediction of stall on-set. 


Although very sophisticated, the notable weakness of all 
viscous codes is that they must rely on empirical models to 
describe transition and turbulence. Quite complex  turbu- 
lence models have been developed in the past twenty years. 
Bradshaw (Ref. 18), Launder and Spalding (Ref. 19), Rotta 
(Ref. 20), Mellor and Herring (Ref. 21), Launder, Reece, and 
Rodi (Ref. 22) all introduced complex methods to solve tur- 
bulent shear stresses using sets of partial differential 
equations. Although superior to Prandtl's mixing length 
theory in their generality, none of these methods results in 
an adequate quantitative description of the turbulent shear 
layers in cascades. Shear layer transition has also been 
widely investigated. Prediction methods for boundary layer 
transitions were developed by Dunham (Ref. 3), McDonald and 
Fish (Ref. 23), and Forest (Ref. 24). Abu-Ghannam and Shaw 
(Ref. 25) measured the transitional boundary layer in a wind 
tunnel to show the effect of free-stream turbulence, pres- 
sure gradient, and flow history. Blair (Ref. 26) measured 
the effect of free-stream turbulence in favorable pressure 
gradients as did Wang, Simon, and Buddhavarapu (Ref. 27). 

Although the measurements show good agreement with 
transition onset predictions, transitional length predic- 
tions lack good quantitative agreement. Prediction methods 
for free shear layer transition, such as occurs over laminar 
bubbles, vere developed by Horton (Ref. 28) and Roberts 


(Ref. 4). Measurements of the separation bubble taken by 
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Walker (Ref. 8) broadly support Horton"s semi-empirical 
model although some aspects of his theory require signifi- 
Gane modifications. A review of the nature of transitional 
processes in turbomachine blading and existing methods for 
prediction of the extent of transition was given recently by 
Walker (Ref. 29). 

In summary, inviscid numerical codes work well. Design 
point conditions in a cascade can be analyzed successfully 
because the boundary layers are relatively thin and remain 
attached, and adequate techniques have been developed to 
model short laminar bubbles and blade wakes. Off-design 
conditions however, are characterized by thickening bound- 
ary layers and, at some threshold, detachment moving forward 
from the trailing edge. The aerodynamic forces and moments, 
the exiting average flow angle and the losses, will then 
depend critically on the thickened boundary layer profile 
and the location of separation. These, in turn, will depend 
critically on the correctness of the description used for 
transition and turbulent growth. Since the physics of 
transition and turbulence cannot yet be reproduced from 
fundamental laws using numerical methods, information to 
complete a numerical simulation must be derived from 
experiments. Such was the goal of the present 


investigation. 
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D. PREVIOUS EXPERIMENTS 

In this section some of the earlier experiments that 
served to influence the present work are briefly reviewed. 

NASA SP-36 (Ref. 30) is an early but thorough report on 
the design of axial compressors. It describes clearly how 
cascade concepts are used in a compressor "design system." 
It contains analysis and correlations of the experimentally 
determined performance of two-dimensional cascades of NACA- 
65 series, C-series, and DCA bladings. Based on the experi- 
mental data, methods to calculate the required blade camber 
angle, design incidence angle, pressure loss, deviation 
angle and, consequently, the blade setting angle, to produce 
specified inlet and outlet velocity vectors, are described. 
Information gained from many, systematic cascade wind tunnel 
experiments was reduced to correlating equations and data 
for blading design. As a result, an axial compressor can be 
designed successfully using the information given in NASA 
SP-36. What is not immediately apparent, however, is that 
all the information qiven refers strictly to che desiqn 
point and can be used to predict effects of only very small 
incidence changes from design point incidence. Also, the 
great care necessary in the conduct of cascade tests is not 
evident in SP-36 since the report deals only with the corre- 
lation of cascade results. The many difficulties involved in 
obtaining inlet flow uniformity, blade passage periodicity 


and span-vise tvo-dimensionality vere encountered and 


J 2 


overcome by the early NASA investigators (e.c., Ref. 31). 
27527 data from strictly two-dimensional flow, requiring 
27 11 suction to produce an "axial-velocity-density- 
ratio" (AVDR), equal to unity, vere used in the design 
correlations. 

In a compressor, AVDR # 1 implies stream surface con- 
traction as it passes through the blading. In a cascade 
wind tunnel, AVDR # 1 for the center-span stream surface is 
the result of boundary layers and corner flow on the side- 
walls at the ends of the blades. Starken, Breugelmans and 
Schimming (Ref. 32) investigated the effect of ADVR on cas- 
cade parameters such as pressure coefficient distribution, 
exit flow angle, loss coefficient and diffusion factor, by 
varying boundary layer suction at different inlet angles. 
They showed the importance of AVDR corrections when compar- 
ing prediction with experiment in two-dimensional cascades. 

Following the early studies of viscous effects leading 
to compressor design techniques (Ref. 30), relatively few 
detailed studies of boundary layers on airfoils in cascades 
have been reported. Evans (Ref. 33) compared measurements 
of suction surface boundary layers on cascade airfoils with 
Similar measurements made at the mid-chord of a compressor 
Stator blade using hot-wire instrumentation. He showed that 
the boundary layer growth on the stator blade was signifi- 
cantly greater than on the cascade blade. Evans explained 


the phenomenon by invoking the unsteady nature of the 
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oundary layer development on the stator blade. This 
information is relevant in applying cascade data to compres- 
sor design and in making 2-D CFD codes applicable to the 
compressor environment. Shama et al (Ref. 34) made hot-film 
and hot-wire measurements in a wind tunnel simulating the 
two-dimensional boundary layer on the suction surface of a 
turbine airfoil, in order to validate computational methods 
and Ə. models. They showed that the McDonald-Fish 
turbulence model [Ref. 23] gave reasonable predictions in 
attached transitional boundary layers. Walker (Ref. 35) 
made hot-wire measurements of the turbulent boundary layer 
on the suction surface of an axial compressor stator blade. 
Walker compared the experimental results with then-current 
calculational methods to assess their reliability. More 
recently Deutsch and Zierke (Ref. 36) used a one-dimensional 
LDV system to measure the boundary layers on the suction and 
pressure surfaces of a DCA blade in a cascade at 5 degrees 
incidence angle. Their work served as a useful reference 
for the present study since it included various other 
measurements such as surface pressures, flow visualization 
to determine on-set of transition and separation, near wake 
measurements using the LDV, and far wake measurements usina 
a five-hole probe. 
Several investigators have reported wake measurements 
made in cascades. Lakshminarayana and Davino (Ref. 37) 


measured stator blade wakes in an axial flow compressor. 
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They drev conclusions concerning Sinto ty decay 
characteristics, wake width, and pressure distribution in 
the trailing edge region, near and far wake. They observed 
substantial enough differences between compressor stator 
wakes and cascade wakes to warrant the development of a 
three-dimensional analysis. Ravindranath and 
Lakshminarayana (Ref. 38) reported similar measurements made 
in the wakes of a compressor rotor blade. Hobbs et al (Ref. 
39) examined the wakes of CD blades in a cascade with the 
intent of improving the modelling of wake effects then using 
inviscid analysis. The experiments included measurements of 
the near and far wakes and of the trailing edge boundary 
layers using hot-film, Kiel, and five-hole pneumatic probes. 
An outgrowth of the latter study was an experiment by 
McCormick and Paterson (Ref. 40) to simulate the CD blade 
trailing edge and wake flows on a large scale. A single one 
inch plate vas used in a model wind tunnel, the walls of 
which were adjusted to simulate mid-passage streamlines. 
Qualitative similarities were found in the wake behavior 
measured in the present study and McCormick and Peterson's 
detailed observations for a single, simulated test condition 
mer a specific, but different, CD blade design. In 
particular, the strong asymmetry in the near-vake velocity 
profiles was observed in both studies. This asymmetry was 
previously reported by Dreon (Ref. 41) based on pneumatic, 


five-hole cone probe measurements of the CD blade wakes in 
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the present cascade, The present vork includes a comparison 
of the results of LDV measurements vith Dreon"s calibrated 
pneumatic probe data, 

The program at the Turbopropulsion Laboratory of the 
Naval Postgraduate School leading to the present study in- 
volved experimental measurements in the lov speed cascade 
vind tunnel first vith DCA blading, and then vith CD blading 
vhich vas designed to replace the DCA section at the mid- 
span section of the stator of a NASA transonic stage, 
Following several preliminary studies, Himes (Ref. 42) re- 
ported the blade element performance of the reference DCA 
cascade design. McGuire (Ref. 43) used flow visualization 
techniques to examine the transition and boundary layer 
separation on the DCA blades. Koyuncu (Ref. 44) replaced 
the DCA blades with the "equivalent" CD blades and evaluated 
their performance over a wide incidence range. Sanger and 
Shreeve (Ref. 45) compared the CD cascade results with 
numerical design and analysis code predictions. The two 
series of blading tests included surface pressure measure- 
ments, calibrated five-hole probe inlet and wake measure- 
ments, and surface flow visualization using tufts, oil smear 
and china clay techniques. A lover loss coefficient and 
wider range were measured with the CD than with the DCA 
blading. Measurements compared reasonably well with results 
of non-interactive numerical codes at the design incidence, 


except for the pressure distribution and the identification 
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of separation near the trailing edge, The prediction of 
losses at off-design incidence angles was extremely poor, 
and the viscous calculations were found to require the input 
of boundary layer information which was not resolved 
adequately in the experiments. 

Thus, the present study capitalized on the previous 
experimental knowledge, but was focused on providing the 
precise viscous flow field information necessary for the 
calibration of strong-interaction or fully viscous analysis 
methods required in the prediction of off-design performance 


and compressor blade stall. 
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III. TEST FACILITY”AND Pİ LEUNELTƏECİL 


A. RECTILINEAR CASCADE 

A schematic diagram of the rectilinear subsonic cascade 
wind tunnel facility is shown in Figure 2. A more detailed 
description of its design and operation was given by Rose 
and Guttormson (Ref. 46). Flow inlet conditions were inves- 
tigated in detail by McGuire (Ref. 43). While uniform on 
average, the inlet flow (far upstream) contains periodic 
wakes due to variable inlet guide vanes which are spaced one 
inch apart in the blade-to-blade direction. The test sec- 
tion downstream of the guide vanes is shown in Figure 3. 
The test section dimensions are slightly different from 
those given by Koyuncu (Ref. 44) and Dreon (Ref. 41). A 
broken line in Figure 3 represents the optical window used 
in the present study. The original 7/8" thick plexiglas 
window was replaced by a 1/4" thick plexiglas window to 


reduce light transmission energy losses. 


B. CD BLADING AND CASCADE CONFIGURATION 

The CD test blades and cascade geometry were scaled from 
the design of the mid-span section of a CD compressor stator 
blade row (Ref. 5 Vand Ref. 45). Twenty aluminum blades 
were machined with a span of ten inches to fit the 60" XN 


test section of the rectilinear cascade. A single blade is 
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shovn in Figure 4. The coordinates of the profile are given 
in Table 1. The qeometrical parameters for the cascade are 
given in Table 2. A single test passage, bounded by two 
adjacent blades as seen through the optical window is shown 
M Figure 5. The coordinates of the test passage are given 
in Table 3 in the coordinate system of Figure 5. 

To reduce back reflections from the blade surface four 
blades were specially treated. Two blades were polished, 
and the other two were black anodized. The latter allowed 
flow measurements in the boundary layer closer to the blade 
surface. 

Pressure taps were provided in two other blades. The 
coordinates of the pressure taps are shown in Figure 6. The 
locations of the different blades in the cascade are shown 
Bü Tigure 3. At each setting of the inlet air angle, the 
blade pair defining the test passage shown in Figure 5 were 


the 7th and 8th blades from the left in Figure 3. 


eR INSTRUMENTATION 
er Doppler Velocimeter (LDV) 

The two-dimensional LDV measurement system consisted 
of four major subsystems: (a) the laser and optics, (b) the 
data acquisition system, (c) the traverse mechanism, and (d) 
the seeding probe. A photograph of the LDV equipment is 


shown in Figure 7. 
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İNCH 
Controlled Diffusion Blade Profile. 


Fıqure 4. 


TABLE 1, | 


PEST TABLADE COORDIYATES (INCHES) | 


X Y (PRESSURE SIDE) Y (SUCTLOM SIDE) 
250 0.045 0.045 
0.022 ae 0.084 | 
15057 0.002 === | 
2222 0.044 0.196 | 
0.444 oz 01 0.307 | 
0.666 0.155 0.403 | 
0.888 02209 0.488 | 
1:110 0.255 0.561 
il 332 0.299 üə 
1.554 0.330 0.663 | 
55776 0.350 0.691 | 
1.998 0.359 0.705 | 
0220 0.359 0.708 
2.442 0.352 05701 
2.664 0.342 0.681 
2.886 ORT 0.650 
3.108 О 0.610 
25330 0.301 0.563 
3.552 0.281 0.510 
27774 üə 0.453 
3.996 ün 0.393 
4.218 ünü Peco 
4.440 0.146 02270 
4.662 0.089 0.208 
4.884 0.019 0.145 
4.925 i — –—- 
5:7 y DEM | 
5.010 0.062 0.062 
(Coordinate System as in Figure 4) 
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| TABLE 2 | 


GEOMETRICAL PARAMETERS OF THE CASCADE 


Number of Blades 20 

Blade Spacing 5.0: 

Chore 5077 

Solidity 1:07 

Trailing Edge Radius 0/062 

Setting Angle 14.2 + 0.1 deg. 
Staqqer Anqle 14.4 £ 0.1 deq, 
Span 10207 





NOMINAL TEST CONDITIONS 





m 20 oR 

Eti 1,032 Atmos. 
P2 1.00 Atmos. 
My ün 

Rec 700,000 
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| 
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: | 
ickness 7* Chord | 

| 
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| 

| 
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i 1 2 3 4 5 6 7 8 
у(") -6.292 -4.792 -4.542 -4.292 -A.042 -3.792 -3.292 -2.792 
9 10 11 12 15 14 15 16 17 18 19 
ил 192 1.292 =0.792 -0:5 -0.25 0.0 0.262 0.362 0.678 1.062 


e 





Y (INCH) 
-2 


1.9 2 26 3 3.5 4 


0.6 1 
Х (INCH) 


Figure 5. Two Adjacent Blades and Test Passage Geometry, 
Showing Survey Stations for LDV Measurements. 
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TABLE 3 


COORDINATES OF THE TEST PASSAGE BETWEEN ADJACENT BLADES 


VERTICAL HORIZONTAL HORIZONTAL 
M -4.7905 T 6 1.7614 
ae -4.7806 -122594 1.7426 
3. -14:7730 -1 1.7364 
4. 247077 71.2674 1479320 
5 -4.7619 i NS J. 73501 
6. -4 7596 =1.2/90; 137 
7. -4.7509 ~ Vee? 25 17277 
8. -4.7407 212721 J. 7279 
9 24273559 EVI qo S 
10. -4.7310 21226995 1.7304 
Ji. -4.7245 -1.2666 1.7334 
12, =—1.57151 –1.1848 1.8152 
Фу -4.3741 =1. 0747 1:9253 
14, -4. 1725 05775 25055 
15. -3 0 02062.) 29013977 
16. =. 70668 -0.7610 222590 
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47. 0.0007 0.0621 3.0621 
48. EDS 0.0646 3.0646 
ES. UIN 0.0708 3.0708 
50. =0 54 36 0.0761 3.0761 
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COORDINATES OF THE TEST PASSAGE BETWEEN ADJACENT BLADES 
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(CONTINUED) 


HORIZONTAL 
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a. Laser and Optics 


A TSI Model 9100-7 (backscatter) LDV system was 


used. A top view of the system is shown in Figure 8. The 
laser was a Laxell five-watt argon-ion laser. A prism 
separated the laser output beam into different colors. Two 


colors, green (514.5 nm) and blue (488 nm) were selected by 
the optical components. The two beams were centered and 
split into a four beam arrangement to measure two velocity 
components at right angles to each other. Two Bragg cells 
shifted the frequency of one beam in each pair to allow 
measurements of reversed flows. The four beams then passed 
through a divergence section which improved the dimensions 
of the measuring volume. The end lens produced a focal 
length of 762 mm for the four beams. Two photo-detectors 
collected the scattered light after it passed through 
the same optics. Table 4 contains a summary of the charac- 
teristics of the LDV system. 
b. Data Acquisition 
The data acquisition system consisted of two TSI 
Model 1990 counter-type signal processors in which the siq- 
nals from the photo-detectors were processed and transformed 
into voltages proportional to the Doppler frequencies. An 
interface unit mounted on one of the counters transferred 
the information to a Hewlett-Packard 1000 series computer 


that controlled the data acquisition and performed 
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Figure 8. LDV System. 
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TABLE 4 


LDV CHARACTERISTICS 


WAVELENGTH 

FRINGE SPACING 

FOCAL LENGTH 

NUMBER OF FRINGES 

HALF ANGLE 

MEASURING VOLUME DIAMETER 
MEASURING VOLUME LENGTH 


5 EOUENCY SHIFT 


SIGNAL PROCESSOR 


COUNTER SETTINGS 


| 
| 


Blue - 488 nm, Green - 514.5 nm) 


Blue - 4.51 um, Green - 4.76 um| 


762 mm 
28 
5.10” 
133 um 
25) mM 


zkEuz3250740 HUz 1201,2,5, 
sequence to 10 MHz 


Counter Type (TSI Model 1990) 
Filters - As Required 

Gain - As Required 

Single Measurement per Burst 
8 Cycles per Burst 
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preliminary analysis of the data. The software "DRP3" was 
provided by TSE: 
c. Traverse Mechanism 

The laser and optics, mounted on an optical 
bench, were installed on the bed of a milling machine modi- 
fied to serve as a traversing mechanism. The bed was manu- 
ally traversed in x,y, and z directions, and the x and z 
coordinates were recorded manually using digital-electronic 
readouts. 

d. Seeding 

The flow was seeded with commercial olive-oil 
droplets. The seeding system consisted of the particle 
generator and the seeding probe. The particle generator is 
shown in Figure 9. The geometry and dimensions of the seed- 
ing probe are shown in Figure 10. The probe was installed 
upstream of the guide vanes to seed the required streamlines 
such that there was minimum disturbance to the flow and 
minimum oil contamination on the optical window. The 
particle size distribution was measured in preliminary tests 
as a function of the feed air pressure as described in 
Apnendix A. An air pressure of 40 psi was selected to give 
an average particle size of 0.9 microns with a standard 
deviation of 0.45 microns: Ninety-five percent of the 
particles were smaller than 1.8 microns in diameter. These 


characteristics were considered to be satisfactory by the 
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criteria described by Durst, Melling, and Whitelaw [Ref. 
577 OD. 275-285, or Dring (Ref. 48). 
e. Aligning Tool 
An aligning tool was designed to position the 
LDV system measuring volume with respect to the blades in 
the test section. Figure 11 shows the aligning tool in- 
stalled betveen tvo blades. Vhen the measurinq volume vas 
positioned at one of the 0.013" diameter holes the four 
beams vere clearly visible on the far side-vall. 
E Survey Probe 

A United Sensor Model 5711 cylindrical pneumatic 
five-hole probe vas used to measure the inlet conditions 
(velocity and flow angle) to the test section. the probe 
was located 1.8 chord lengths (1.8c) axially upstream of the 
blade row. A Hewlett-Packard Data Acquisition System (HP- 
3032) and an HP Interface Bus (HP-98034 HP-IB) were used to 
collect data from pressure probes and surface pressure taps, 
through two 48-port Scanivalves. The system was controlled 
by an HP-9845A computer. 

3. Cascade Wind Tunnel Operating Instrumentation 

Tunnel (plenum) stagnation pressure was measured 
using a tube suspended in the plenum chamber below the wind 
tunnel. An  iron-constantan thermocouple, similarly sus- 
pended in the plenum, measured stagnation temperature. Two 
rows of static taps (Fig. 3) were connected to water 


manometers, and used to monitor the static pressure 
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27 “ribution into and out of the test cascade, At each in- 
let air  angle,the inlet guide vanes and outlet tailboards 
were adjusted to produce nearly uniform static pressure in 
the blade-to-blade direction both upstream and downstream of 


the test blading. 


D, MEASUREMENT UNCERTAINTY 

Table 5 summarizes the measurement uncertainties. The 
uncertainties of the LDV parameters (K, L, A) were supplied 
by TSI. Fringe spacing was calculated using (Ref. 47) 


A 


15 z sink 


The particle velocity uncertainty vas calculated using 


AV = ne + At = 7: + VAt 
V Af E de вар 
where At = clock accuracy = 1 n-sec. 
t = time to cross 8 fringes = 8df£/V 


Since the particles do not follow the flow precisely, 
the flow velocity uncertainty is always than the particle 


velocity uncertainty. 
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(for V=10 m/sec.) 
(for V=100 m/sec.) 






UNCERTAINTY 


20 arcsec. 


ГУ ЕТ РЕР ТМЕМТАГҺ РБОСЕООЈЕЕ 


A, OVERVIEW 

First, in order to exercise and verify all aspects of 
the LDV measurement system, an experiment vas conducted in 
the cascade vind tunnel to measure the boundary layer on a 
flat plate and to compare the results vith analytical solu- 
tions and other published data. A detailed description of 
the flat plate experiment is given in Appendix B. The 
velocity measurements which were obtained at low Reynolds 
number showed very good agreement with the Blasius solution 
for a laminar boundary layer on a flat plate. The turbu- 
lence distribution measurements also showed good agreement 
with other published data. 

The CD blading was then installed, as shown in Figure 3 
and detailed flow measurements were made for three different 


inlet air angles, namely; 


ıl, a 1: - the design condition, 
2. eee 5:50” — a moderate positive incidence angle, and 
7 :) = 462 = ta high incidence angle, near to stall. 


At each inlet air angle measurements were made covering 
four regions of the flow field, namely; 
1. the inlet flow field 0.3 chords upstream of the blade 


leading edge, 
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2. the passage between the 7th and 8th blades, to map the 
inviscid field, 

3. the 7th blade pressure and suction surface boundary 
layers, to map their growth, and 


4. the wake downstream of the 7th blade. 


By SET-UP PROCEDURES 
1. Wind Tunnel Adjustments 
At each inlet anale, special attention was given to 
setting che tunnel to obtain a uniform and correct inlet 
flow angle. The inlet side-walls were first set to the re- 
quired angle. Then the quide vanes were set approximately 
based on previous experience. The flow angle distribution 
at the inlet was measured with a five-hole probe and then, 
for the test passage only, with the LDV. The guide  vanes 
were then adjusted to give the correct flow angle distribu- 
tion guided by the probe and LDV measurements. Static pres- 
sure distributions were observed using water  manometers to 
verify uniformity of pressures at the inlet and outlet of 
the test section. 
(a) Inlet Flow Quality 
For conditions similar to those in the present 
study, McGuire (Ref. 43) and Dreon (Ref. 41) made probe 
measurements of the inlet velocity field 1.8 chords upstream 
of the test blades and obtained acceptably uniform distri- 


butions of the velocity and the flow angle. In the 
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present study the LDV system vas used to measure the inlet 
conditions to the test passage (between the 7th and 8th 
badges) 0.3 chords upstream of the test blading. (This was 
as far upstream as the window would allow.) The results, 
plotted in Figure 12, showed nearly uniform velocity and 
flow angle at each of the three inlet air angles. The 
deviation from total uniformity is likely to be the upstream 
influence of the blades. 
(b) Periodicity 
At similar conditions to those in the present 
study, Koyuncu (Ref. 44) and Dreon (Ref. 41) made probe 
measurements of the far wake to verify periodic behavior of 
the cascade.  LDV measurements presented in Figure 13 showed 
Similar results for the 7th and 8th blade wakes. 
7  LDV Optical System Alignment 
Periodic LDV alignments were carried out regularly 
to ensure proper beam crossing and fringe patterns. Special 
optical configurations were set for different parts of the 
experiment. Usually it was preferred to orient the beams 
such that the velocity vector lay at about 45 degrees to the 
blue and green beams. Thus, as shown in Figure 14, when 
testing in the lower section of the cascade (Stations 1-5 in 
Figure 5) the beams were arranged in a vertical-horizontal 
arrangement, whereas in the upper part of the cascade 
(stations 6-19 in Figure 5) the beams were rotated by 45 


degrees. Most measurements were taken with the Bragg cells 
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installed, applying proper frequency shifting. However, 
where no reversed flows could be present, frequency shif:ing 
was not needed. The inlet flow field (reference velocity 
field) and the lower section of the passage were measured 
without the Bragg cells installed. 
3. Traverse and Measuring Volume Alignment 

A rectangular coordinate system was chosen where 
the origin lay on the center-of-radius of the 7th blade 
trailing-edge. Figure 5 shows the test passage and the 
coordinate system. Aligning the measuring volume in the 
coordinate system for the passage flow measurements was 
carried out using the special tool shown in Figure ll. The 
table was adjusted such that the four beams crossed at a 
specific hole, and the image was clearly visible on the rear 
wall of the tunnel. The hole coordinates were then set on 
the traverse mechanism digital electronic readouts. The 
measuring volume was aligned for wake measuremerts in the 
same way. For boundary layer measurements the procedure was 
different. First. the optical bench was yawed on the mill- 
ing table, 3.5 degrees with respect to the span-wise 
direction. Then the measuring volume was positioned at the 
proper level (y axis) using the aligning tool. Finally, the 
measuring volume was traversed horizontally (x direction) 
until the four beams crossed on the blade surface. The 
crossing was observed with the naked eye and an accuracy of 


0.002" in the adjustment was achieved with ease. 
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577 LDV SURVEY MEASUREMENTS 
1. Free Stream Measurements 
Free stream measurements were taken 0.3 chords up- 
Stream of the blade row leading edge (station 1 in Figure 
5). The velocity and flow angle distribution (in Figure 12) 
were uniform within + 1% and + 0.75° respectively, despite 
the proximity of the survey station to the blades. These 
measurements provided the reference velocity and reference 
turbulence level for subsequent surveys at the same test 
condition, as described in Appendix C. 
2. Passage Measurements 
Passage measurements were taken at 14 levels through 
the passage (stations 2-15 in Figure 5) to measure the 
inviscid behavior of the flow field. 
3. Boundary Layer Measurements 
Boundary layer measurements were taken at several 
stations on the suction and pressure surfaces. The stations 
were coincident with those of the passage measurements. The 
first survey was made at the trailing-edge center-of-radius. 
The surveys were then moved progressively upstream towards 
the leading-edge until the boundary layers were too thin to 
be measured accurately with the LDV. The procedure for sur- 
veying the boundary layer was first to position the measur- 
ing volume on the blade surface at the required level. The 
LDV was then traversed normal to the surface (using both x 
and y traversing mechanisms) and data were acquired at 
predetermined intervals. 
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4, Wake Measurements 

Wake measurements were made at four different levels 
downstream of the 7th blade trailing-edge (stations 16-19 
in@Pigure 5). 

5. Two-Dimensionality Checks 

Special attention was given to the question of 
whether the flow field at mid-span was two-dimensional. 
Verification checks included the following: 

a. Surveying the wake at y = 0.678" (0.123c) at tvo 
different spanvise locations and comparing the 
distributions vith the mid-span distribution. 

b. Surveying the suction surface boundary layer at 
tvo different spanvise locations and comparing the 
distributions with the midspan distribucion. 

c. Surveying the pressure surface boundary layer in the 
same manner as the suction surface boundary layer. 

The results of these measurements at each inlet air 
angle are plotted in Figures 15-17. It was concluded 
that the flow was two-dimensional at mid-span to within 
acceptable error. 

6. Non-Dimensionalizinq the Data 

The LDV data vere non-dimensionalized by dividinq 
the reference velocity at the time of  recording. The 
reference velocity was calculated using the plenum pressure 


and temperature as described in Appendix C. The plenum 
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pressure and temperature vere constantly monitored and 
recorded throughout the experiment, 
7. Reynolds Number 
The chord length Reynolds number of the experiment 
was nearly constant (Re, = 740000 + 40000). The variation 
were the result of variations in the plenum pressure and 
day-to-day variations in atmospheric conditions. 
8, Data Reduction 
The data vere acquired using the HP-1000 series 
computer, stored on floppy disks and thereby transferred to 
an IBM 3033 (main frame) for analysis. Three Fortran pro- 
grams were written to analyze the data. A listing of the 
software can be found in Appendix D. The programs include; 


i. CASCADE FORTRAN - Analyzes the data from passage 
measurements. 


ii. BLAYER FORTRAN - Analyzes the data from boundary 
layer measurements. 


iii, WAKE FORTRAN - Analyzes the data from wake 
measurements. 

D. SURFACE PRESSURE MEASUREMENTS 

Surface static pressure measurements were taken using 
the partially-and fully-instrumented blades. The partially- 
instrumented blade was used to measure the trailing- E 
static pressure. The fully-instrumented blade was used to 
measure the pressure distribution over the blade at the 
highest incidence angle for which previous data were not 


available. 
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Vie” RESULTS AND DISCUSSION 


In this chapter the measured data are presented, and 
discussed, dealing with the different regions in the flow 


field one at a time. 


A. 81 FLOW FIELD MEASUREMENTS 

The following parameters were measured for each inlet 
air angle; two velocity components, u, v, turbulence level 
€, plenum pressure and temperature Pto, Tto and atmospheric 
pressure Patm: The velocity measurements throughout the 
flow field were subsequently normalized with respect to the 
inlet reference velocity at corresponding inlet air angle. 
The inlet reference velocity for each individual LDV mea- 
surement was derived from the inlet flow field measurements 
and local plenum conditions as shown in Appendix C. This 
procedure removed the effects of varying supply conditions 
while surveys were conducted, and gave results in a form 
suitable for comparison with code predictions. 

Inlet mach number was about 0.25. Inlet velocity was 
about 85 m/sec. 

The free stream turbulence level was measured to be 
e= 1.4+ 0.2% at each test condition. (Free stream turbu- 
lence level was calculated using only u' and v' as defined 


in the list of symbols.) 
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Be FLOW FIELD IN THE PASSAGE 

Streamlines of the flow in the passage for the three in- 
let angles are shown in Figures 18-20. The streamlines were 
generated in the following manner.  Volumetric flow was cal- 
culated for each level by spline interpolation and integra- 
tion of the vertical component of the velocity (v). The 
field was marked every 1/6 of the volumetric flow (to 
generate five streamlines). The flow angle was drawn at 
each mark using the local velocity components (u, v) 
obtained by interpolation. The flow field was seen to be 
smooth throughout the range of angles. While the streamline 
closest to the blade suction surface moved progressively 
away from the surface as the inlet angle increased, there 
was no obvious indication of significant separation. 

1. VelocitviDistriburrem 

Normalized velocity distributions in the passage for 

three inlet angles are shown in Figures 21-23. A comparison 
between the passage inlet (leading edge station) velocity 
distributions at three inlet angles is shown in Figure 24. 
A comparison of the outlet velocity distributions is shown 
in Figure 25. The distributions of the edge velocity (de- 
fined as the velocity where it reaches its maximum value or 
where turbulence level reaches free stream values) for the 
three angles are shown in Figures 26 and 27. There seemed 
to be similar trends at the three inlet angles, namely; 


a. a velocity peak at the suction surface leading edge, 
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Figure 24. Passage Inlet (Leading Edge) Velocity 
Distribution at Three Inlet Air Angles. 
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Figure 26. Edge Velocity Distribution 
Over the Blade. 
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BETA = 46 DEG. 








1.2 1.4 1.6 1.8 2 


EE 


Us/VREF 


0 0.20.4 0.8 0.8 1 





0 10 20 30 40 50 60 70 80 90 100 
7; CHORD 


Figure 26. (cont.) Edge Velocity Distribution 
Over the Blade. 
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b. 


a significant velocity difference between the suction 
side and the pressure side near the leading edges 
(passage inlet) which decreases downstream and becomes 
almost flat over the aft 15% of the passage, 
an almost flat velocity distribution over the pressure 
surface and a slight acceleration towards the trailing 
edge. 

Other observations were as follows: 
The maximum velocity at the suction side inlet in- 
creased as the inlet flow angle increased. Downstream 
of 15% chord the edge velocity on the suction surface 
decreased as inlet flow angle increased. 
The magnitude of the velocity over the pressure side 
decreased as the inlet flow angle increased, as did 
the nearly common magnitude of the suction and pres- 
sure side edge velocities at the trailing edges. The 
distribution of the edge velocity over the pressure 
side surface varied very little as the inlet angle was 
changed. 
At the two largest inlet angles (43.4°, 46°) the first 
two stations (2 and 3 in Figure 5) on the suction sur- 
face gave a nearly constant edge velocity. This 
flatness in the edge velocity distribution was not 
consistent with the measured surface pressure 
distribution as is discussed later. It is suspected 


that the seed particles could not follow’ the high 


21 


acceleration around the blade leading edge, creating a 
velocity laq. Thus, the leading edge data points on 
the suction surface for the higher inlet angles were 
plotted in Figure 26, but marked out. 

d. At the outlet the boundary layer vas verv thick com- 
pared to the passage width. For the high inlet angle 
(6) = 46°) the boundary layer thickness was more than 
20% of the passage width and 12% of the chord. 

2. Flow Angle Distribution 

Flow angle distribution in the passage for three in- 
let angles are shown in Figures 28-30. A comparison between 
the flow angle distributions across the passage at the blade 
row leading edge (Station 2 in Figure 5) for three inlet 
angles, is shown in Figure 3l. A comparison of the passage 
outlet (Station 15 in Figure 5) flow angle distributions is 
shown in Figure 32. The following were concluded from the 
figures: 

a. Most of the flow turning took place in the forvard 
half of the cascade. 

b. The outlet flow angle distribution was almost identi- 
cal at all three inlet angles. Deviation angle did 
not vary Significantly with inlet angle throughout the 
inlet angle range. 

3. Surface Pressure Distribution 
Pressure coefficient distributions over the airfoil 


for the three inlet angles are shown in Figure 33. The 
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Figure 29-en Plow Angle Distribution in the 
Passage; B1 z 43.49, Re = 740000. 
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Figure 29. (cont.) Flow Angle Distribution in the 
Passage; 8, € 43.45, Re - 740000. 
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Figure 30. Flow Angle Distribution in the Passage; 


B1 z 46?, Re - 730000. 
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Figure 32. 
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Figure 33. Pressure Coefficient Distribution 
Over the Blade. 
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BETA = 46 DEG. 
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Figure 33. (cont.) Pressure Coefficient Distribution 
Over the Blade. 
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distributions follov the same general trends as the edge 


25 city distributions in Figure 26, except at the suction 
surface for the higher inlet angles (43.4°, 46°). If the 
boundary layer is thin the edge velocity and surface pres- 
sure are related using the inviscid, incompressible 


expression 
5 
A comparison betveen the edge velocity calculated using this 
expression and the measured edge velocity is shovn in Figure 
34, Since, as suggested earlier, the seed particles cannot 
handle the high acceleration of the flov at the suction sur- 
face leading edge, the leading edge velocity data points 
were omitted as in Fiqure 26. The slight acceleration of 
the flow on the pressure surface towards the trailing edge 
shows as a slight decrease of the pressure coefficient, to 
Satisfy the Kutta condition within the base flow. 
4. Axial Velocity Density Ratio (AVDR) 

The "referred" passage volumetric flow rate, calcu- 
lated at each station by numerically integrating the distri- 
bution of v/Vref from one side of the passage to the other, 
is shown in Figure 35. Eighteen data points are shown 
plotted corresponding to the levels in Figure 5 but exclud- 
ing the near wake, station 16 at y = 0.262". The broken 


lines are a least-squares fit to the data, The calculated 
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BETA = 40 DEG. 
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BETA — 46 DEG. 
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Volumetric flov is seen to fluctuate close to the passage 
inlet since the velocities near the suction surface vere not 
easily resolved and the velocity (including flow reversals) 
within the suspected separation bubble could not be measured 
at all. Fluctuations away from the leading edge have no 
obvious trends nor interpretation. The magnitude of the 
uncertainty in the result can be gauged from an examination 
of the AVDR. 

The AVDR at each station vas calculated by normaliz- 
ing the volumetric flow rate to the reference inlet volu- 
metric flow rate and assuming that the flow 15 incompres- 
sible throughout. The results are shown in Figure 36. The 
data away from the leading edge are seen to be within 1% of 
the linear distribution shown within the passage. The AVDR 
of the test passage was 1.025 + 0.002 for the three inlet 


angles. 


(Dum BOUNDARY LAYER MEASUREMENTS 
57 “üction Sürface 
The suction surface velocity distributions obtained 


at the three inlet angles are shovn plotted in Figures 37- 


BS. The velocity is shown normalized with respect to the 
edge velocity. Both velocity components, normal and 
parallel to the wall, are shown. The following is a 


discussion of the results: 
a. Measurements inside the suspected laminar sepation 
bubble were impossible to obtain since seed particles 


Sm 
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Figure 36. AVDR Distribution Along the Passage. 
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seed particles apparently did not penetrate into the 
separated region. The presence of The bubble, 
however, was evident in the measurements. At the de- 
sign condition SUN Figure 37) at Station $5 
negative velocity component normal to the wall 
(velocity towards the wall) indicates that the flow 
was approachinq reattachment. The fact that measure- 
ments could not be achieved close to the wall at that 
station, as was possible at other stations, was taken 
to indicate that the station was sliqhtly upstream of 
the reattachment point. Further dovnstream, at 
Station 4, the tvo component velocity profiles in 
Figure 37 indicate that the station vas slightly dovn- 
stream of the reattachment point. Ас 81 = 43.49, ас 
Station 3 the velocity profiles in Figure 38 indicate 
that this same station vas nov above the bubble. The 
folloving profile, at the same angle, indicates that 
Station 4 vas slightly upstream of reattachment. 
Further downstream, at Station 5, the flow was slight- 
ly downstream of reattachment. At 6, = 46°, shown in 
Figure 39, the measurements exhibited the same quali- 
tative behavior. Station 3 was above the bubble. 
Station 4 was also above the bubble. Station 5 was 
close to reattachment and Station 6 was close but 
downstream of reattachment. The data suqqest that as 


the inlet angle was increased the separation bubble 


96 


bubble moved or extended further downstream. te is 

noted that the component of velocity normal to the 

surface was negligibly small at all stations down- 
stream of the indicated reattachment. 

Distributions of the turbulence level through the 

suction side boundary layer are shown plotted in 

Figures 40-42. The following observations are made: 

(1) High levels of unsteadiness were recorded above 
the bubble. At the design condition (8j = 40°, 
Figure 40) turbulence there was 10-12%. At 
higher angles the unsteadiness progressively 
increased, with e = 12-13% at By = 43.4° and 16- 
173 at 81 = 46%. 

(2) The turbulence level dropped to 8% and lower 
after reattachment. 

(3) The peak in the turbulence level profile moved 
away from the vall as the station moved 
downstream. 

(4) Very close to the wall, where the velocity 
gradient was high, high turbulence levels were 
recorded, but were due to velocity broadening 
(Significant velocity variations in the measuring 
volume at the LDV focal point). Those data 
points were, therefore, shown in the plots with a 


cross. 
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Figure 42. 


From the high levels of turbulence above the bubble 
and the observation that all profiles downstream of 
reattachment were clearly not laminar it was concluded 
that transition occurred as a result of the bubble. 
The boundary layer growth is shown in Figures 43-45. 
Figure 43 shows the growth in the thickness of the 
boundary layer along the suction surface of the blade. 
The boundary layer thickness was defined as the dis- 
tance from the wall where the component of the 
velocity parallel to the wall was 99% of the edge 
velocity. Figure 44 shows the displacement thickness 
distributions and Figure 45 the momentum thickness 
distributions, calculated by spline interpolation and 
integration of the boundary layer data. 

At the trailing edge, as the incidence was 
increased, the suction side boundary layer became very 
thick although no separation was detected. For the 
highest inlet angle (8j » 46?) the displacement thick- 
ness was about 6.5% of the passage width and the 
boundary layer thickness was more than 20% of the 
passage width. 

The distribution of the shape factor along the 
Suction surface is shown in Figure 46, After 
reattachment, the shape factor at first decreased 
until a minimum value was reached, and then increased 


as the turbulent boundary layer thickened. The 
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minimum shape factor vas at approximately 30% chord. 
The value at the minimum was approximately 1.7 for the 
@esign condition (8j = 40%) and approximately 1.5 for 
the higher angles (8) = 43.4°, 46°). 

The similarity in the shape factor distribution at 
the two largest inlet angles but clear departure from 
the results obtained at the design inlet angle implies 
a nonlinear behavior in the boundary layer 
development. This would not be expected simply by 
looking at the results in Figures 44 and 45 for the 
displacement and momentum thicknesses respectively. 
Each of the thicknesses, and the overall thickness 
shown in Figure 43, increased steadily as the inlet 
angle was increased. However, Figure 24 and Figure 27 
show that the velocity distribution near the leading 
edge on the suction side behaved somewhat differently 
at the design angle compared with the two higher 
angles. Also, it can be seen in Figure 33 that no 
strong leading edge suction peak was detected at the 
design angle, but one did occur at the two larger 
angles. Unfortunately, the differences in the early 
development of the boundary layer, particularly the 
presence, structure, and effects of the leading edge 
bubble, could not be resolved in the experiment. 

It is noted that the magnitude of the shape factor 


at the design angle increased to 2.2 at the trailing 
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edge, vithout sign of separation. Fortuitously, the magni- 
tude of the shape factor vas less at the off-design, hiqher 
incidence angles. 

Tables 6-8 summarize the thickness data and their 
uncertainties. The method for calculating the uncertainty 
is described in Appendix E. 

2. Pressure Surface 

Pressure surface boundary layer profiles are shown 
plotted in Figures 47-49. Transition from a typical laminar 
velocity distribution to a typical turbulent distribution 
seen, at each angle, as the survey station moves downstream. 
The pressure surface boundary layers are clearly much thin- 
ner than the suction surface boundary layers, especially at 
the highest inlet angle. 

The growth in the boundary layer thickness on the 
pressure surface is shown in Figure 50. The calculated 
displacement thickness distribution is shown in Figure 51 
and the momentum thickness distribution is shown in Figure 
57 The folloving observations are made vith references to 
Figures 50-52. 

a. The inlet angle had a negligible effect on the mea- 
sured boundary layer thickness. However, the higher 
the inlet angle, the smaller became the displacement 
and momentum thicknesses. The very small changes and 
the trend were due to the fact that the edge velocity 


on the pressure surface changed slightly from one 
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inlet angle to the other, and the higher the inlet 
angle the lower was the edge velocity. (See Section 
5. 

A typical transitional boundary layer on a flat plate 
is characterized by a progressive reduction of the 
displacement thickness until the boundary layer has 
become turbulent (Ref. 49). This trend is seen in 
Figure 51 at the two highest inlet angles (43.4°, 
46°), and the minimum is seen to occur at about 50% 
chord. At the design angle (40°), the displacement 
thickness grew steadily from about 30% chord. 

The momentum thickness, although smaller in magnitude, 
behaved very similarly to the displacement 
thickness. 

All thicknesses decreased slightly over the aft 5% of 
chord. This was the result of the local acceleration 
of the outer flow over the pressure surface near the 
trailing edge, which is clearly seen in Figure 27. 

The shape factor distribution is shown plotted in 
Fiqure 53. Three regions are clearly evident at each 
inlet angle. The forvard region, vhere the shape 
factor decreases sharply, is indicative of a transi- 
tional boundary layer, The central region, with 
a nearly tlat distribution at H= 1.7 -1.8, is 
typical of an attached turbulent boundary layer. The 


aft region where the shape factor is seen to increase 
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quite suddenly indicates that the flov is about to 
separate over the trailing edge curvature. 

Data for the calculated thicknesses and the uncer- 
tainties resulting from the positioning of the probe 
volume at the pressure surface, are summarized in 
Tables 9-11. A further uncertainty is caused by the 
need to interpolate in the integration between the 
measurements and zero velocity at the surface. A 
larger fraction of the layer is involved in the inter- 
polation when the layer is thin, as is the case on the 
pressure side. 

The onset of transition also could not be resolved 
in the pressure side measurements since the laminar 
boundary layers were simply too thin to be measured 
properly by the LDV system. The edge of the transi- 
tional zone, however, could be roughly estimated after 
an examination of the turbulence intensities. The 
boundary layer turbulence level distributions are 
shown plotted in Figures 54-56. As vas the case on 
the suction surface, high instability vas evident in 
the transitional boundary layer region covering at 
least 30% of the chord (e = 103-123). For the design 
eagle (8 = 40°) the boundary layer was fully turbu- 
5 117 :7105— 00 “chord. For the higher angles (8) = 
43.4°, 46°) the boundary layer was fully turbulent 


at 60% chord. Full transition at the higher angles 
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vas delayed due to the fact that the edge velocity 
over the pressure side of the airfoil vas lover than 
at the design angle. 

The fully turbulent boundary layer (beyond 705 
chord) had a lower peak turbulence level (e = 62). 
The distribution of the turbulence level in the fully 
turbulent boundary layer on the pressure side was dif- 
E Nt than that on the suction surface since the peak 


in the turbulence level occurred closer to the wall. 


D. WAKE MEASUREMENTS 

Time Averaged Velocity 

Wake velocity distributions are shown plotted in 
Figures 57-59. The following observations are made with 


reference to the Figures: 


Ae 


The variation of the minimum of the vertical component 
or ehe velocity (which is im the direction of the in- 
tended downstream flow at the design condition) as a 
function of the distance from the trailing edge, is 
seen to be typical of wake decay. At design condition 
5:07”) EOL y = 0.2627 ethoeminimum velocity was 
negative, i.e., the flow was reversed due to separa- 


tion and recirculation in the base of the trailing 


edge curvature. For y = 0.362" the average minimum 
velocity was positive, but small. Bos y= 0.678" the 
minimum velocity was about 28% and for y = 1.062" 


about 40% of the inlet reference velocity. 
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At the higher inlet angle (8; = 43.4°) hc ҝу. 


trends were seen but the magnitude of the minima were 


lower. At the highest inlet angle (3, = 46°) the 
results were somewhat different. At y = 1.062" the 
minimum vas again lover than at 3) = 43.4° and for y = 
0.678" the minimum vas the same, But the near vakes 
at 8) = 46° showed higher minimum velocities ri ma m 
the lower inlet angles. Although, some data samples 


showed negative velocities in the near wake at fj = 
46° the time averaged mean velocity was everywhere 
positive. 


At design inlet angle (8) = 40°) the wake profile was 


almost symmetric. The suction surface side of the 
wake was Slightly thicker than the pressure side. At 
higher inlet angle (8j = 43.4?) the wakes were clearly 


not symmetric, the suction side being much thicker 
than the pressure side. At the highest inlet angle, 
the wake profile had pronounced asymmetry. 

The distribution of the horizontal component of the 
the velocity through the wake is seen to follow a wave 
shape. Traversing from the pressure side towards the 
center of the wake resulted in a gradual increase of 
the horizontal velocity component until a maximum was 
reached. From there, there was a sharp decrease to 
negative values, passing through the center of the 


wake. Traversing further, the horizontal velocity 
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component gradually increased to small positive 
values, This vave shape vas found through all 
measured wakes except at the design angle in the most 
far wake where the profile was somewhat smeared. 

A comparison of the wakes at y = 0.678" at three in- 
let angles is shown in Figure 60. The decrease of the 
minimum of the vertical component as the inlet angle 
increased is clearly evident. The pressure surface 
side of the wake is seen to have changed only slightly 
with inlet angle. In contrast, the suction surface 
side of the wake is seen to have thickened 
significantly as the inlet angle was increased. 

The corresponding distributions at y = 0.678" of the 
flow angle, derived from the two components of 
velocity, are shown in Figure 61. The flow angle 
exhibits the same qualitative behavior as the horizon- 
tal component of the velocity. The pressure surface 
side of the wake did not change significantly with in- 
let angle, whereas the suction surface side was 
changed sianificantly. 

2. Turbulence 
Distributions of the turbulence level measured in 
the wakes are shown plotted in Figures 62-64. The following 
observations are made with reference to the Figures. 
a. The pressure surface side of the wake vas charac- 


terized by higher levels of turbulence. The peak 


143 


0.8 
x 








© i 
o: . 
b 
(xy € 
=“ ә ас иә әта си E ә ә R әти әта ии 
Ҹ | — BETAz40... 
S x BETA=43.4 
x ВЕТА=46 
x 5 
e 1 l 
N 
n 
-1.5 -1 -0.5 0.5 1 1.5 


0 
X (INCH) 


Figure 60. Velocity Distributions at y = 0.678" 
(0.123c from Trasillinoğin lr 


144 


FLOW ANGLE 
0 20 40 


-20 


LEGEND 


o BETA-43.4 
6 BETA=46 





=40 


-1.5 zi -0.5 0.5 i 1.5 


0 
X (INCH) 
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level reached 15% in the near wake (y = 0.262") 
and dropped to 12% at the most downstream station 
(y = 1.062") for all three inlet angles. 

The suction surface side of the wake was characterized 
by lower turbulence levels. The peak levels were 
about 8%-9% at all inlet angles and positions 
downstream of the trailing edge. 

memenewscesign angle (8) ” 40”) turbulence level dis- 
tributions characteristic of the two blade trailing 
edge boundary layers remained defined in the near 
wakes (v = 0.262", 0.362"). They are seen to be 
smeared out at the two downstream stations, and at y = 
1.062" the distribution was smooth and nearly 
symmetric. 

At the higher inlet angles (8) = 43.4°, 46°) the tur- 
bulence level distribution characteristic of the tvo 
boundary layers remains evident at all stations, 
particularly at the highest inlet angle. 

A comparison of the turbulence level distribution 
for the three inlet angles at y = 0.678" is shown in 
Figure 65. The pressure surface side of the wake is 
seen to be almost identical for all three inlet 
angles. On the other hand, the suction surface side 
shows the wake broadening, but with little change in 
the peak turbulence level, as inlet angle was 


increased. 
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Figure 65. Turbulence Level Distributions at y = 0.678" 
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A. 


VI. ANALYSIS OF THE RESULTS 


COMPARISON WITH PREVIOUS EXPERIMENTS 


The following is a summary of the main parameters and 


physical features identified in the previous chapter in 


comparison with previous experience. 


1, 


25 


The free stream turbulence level was 1.4 + 0.2%. 

The AVDR did not change significantly with inlet angle 
variation. An AVDR of 1.025 was calculated between 
stations 0.3 chords axially upstream and 0.24 chords 
downstream of the blade. This was consistent with 
Dreon (Ref. 41), who calculated AVDR using measure- 
ments from a five-hole probe to be 1.05 + 0.01 between 
stations two chords axially upstream and 1.65 chords 
downstream of the blade. 

Measurements of the edge velocity agreed quite well 
with the velocity calculated from measurements of the 
surface pressure in those regions of the flow where 
the surface curvature was moderate and the boundary 
layer was relatively thin. The edge velocity and cal- 
culated (inviscid) surface velocity would be expected 
to be different where the surface curvature is signif- 
icant and the boundary layer not sufficiently thin. 

A separation bubble was identified on the suction 


surface by the LDV measurements. A rough estimation 
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Or the reattachment point using the LDV data 
was  5$-15$ of chord. Sanger and Shreeve (Ref. 45) 
reported that flow visualization with china clay lo- 
cated reattachment somewhat closer to the leading 
edge, at 4%-7% of chord. Their results showed the 
same trend as did the LDV data, i.e., the reattachment 
point moved downstream as inlet angle was increased. 
Sanger and Shreeve's results were obtained at a chord 
Reynolds number of 340000. The lower Reynolds number 
may explain the quantitative difference. In both the 
present and previous experiments with china clay, it 
was concluded that transition took place in the free 
shear layer, above the bubble. 

The turbulent boundary layer on the suction surface 
vas found to be fully attached throughout the inlet 
angle range. Sanger and Shreeve (Ref. 45) used the 
china clay method to detect separation. Their find- 
ings were different. At 8) ” 39? separation vaza. 
tected at 95% Chord: At 8j = 43° separation was at 
90% chord and at Bi = 46° separation was around 87$ 
chord, Once again, the difference in Reynolds number 
may explain the apparent disagreement. However, an 
alternative explanation can be that the china clay 
technique gave erroneous results in this particular 
region of the flow field. The vertical orientation of 


the suction surface near the trailing edge allowed the 
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liquid in the china clay to creep forward against very 
low dynamic pressures within the suction surface 
boundary layer. the interpretation that the drying 
pattern indicated separation could be in error under 
these circumstances. LDV measurements at lover 
Reynolds number are needed to resolve the question. 
Inlet angle variations had relatively little effect on 
the boundary layer distribution on the pressure 
surface, although increasing the inlet angle (at posi- 
tive incidence anales) delayed transition and caused 
the turbulent boundary layer at the trailing edge to 
be measurably thinner. This is consistent with 
previous experience (Ref. 30). 

A comparison between the LDV wake measurements in the 
present study and measurements obtained by Dreon (Ref. 
41) can be seen in Figures 66-69. Dreon used a cali- 
brated five-hole pneumatic probe to measure the 
velocity and flow angle distributions at six stations 
downstream of the trailing edge. The comparison is 
made with Dreon's data obtained at 0.123 chords down- 
stream of the trailing edge for the design inlet angle 
(40?) and one positive incidence angle (8, = 43.4?). 

At the design angle shown in Figure 66, there was 
Seen to be qood agreement between the two velocity 
distributions. There was a qood match in the location 


of the center of the wake, the minimum velocity, and 
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Figure 66. Wake Velocity Distribution; 
81 = 4077 
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Figure 67. Wake Flow Angle Distribution; 
81 = 407. 
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Figure 68. Wake Velocity Distribution; 
81 = 4305457 
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Figure 69. Wake Flow Angle Distribution; 
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in the velocity throuqh the tvo sides of the vake, 
There vas some difference, hovever, betveen the flov 
angle distributions shown in Figure 67. Larger excur- 
sions of the flow angle through the wake shear layers 
were detected in the LDV measurements than were indi- 
cated in Dreon's probe measurements. At the higher 
inlet angle shown in Figure 68, the two velocity dis- 
tributions agreed well in the location of the center 
of the wake and in the distribution through the pres- 
sure surface side of the wake. There was disagreement 
in the magnitude of the minimum velocity and in the 
distribution through the suction surface side. The 
LDV system measured a lower minimum and a thicker 
wake. The flow angle distribution in Figure 69, 
showed again quite a significant disagreement although 
the general trends were similar. It should be noted 
that Dreon attempted to resolve the uncertainty in his 
measurements of the flow angle in the blade wakes by 
uSing a special probe designed to measure angle (only) 
through the shear layer. Unfortunately, no measure- 
ments were made with the angle probe at the stations 
where LDV data were taken in the present study. Thus, 
the comparisons in Figures 68 and 69 are with data 
from a United Sensor cone probe, with the probe shaft 
quite close to the blade trailing edge. Probe 


interference may explain the differences. 
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8. 


A comparison between some LDV wake measurements in the 


present study and measurements obtained by Baydar 
(Ref. 50) can be seen in Figures 70-71. Baydar used a 
single wire TSI hot-wire system to measure the velo- 
city and turbulence level distributions at three 
stations in the wake. The comparison is made with 
Baydar's data at 0.123 and 0.2 chords downstream of 
the trailing edge for design inlet angle (8, = 40°) 
and the highest incidence angle (8i = 46°). The com- 
parison shows very good agreement between the two sets 
or data. 

An attempt was made to calculate the loss coefficient 
from the LDV data to compare with data reported by 
Sanger and Shreeve (Ref. 45) and Dreon (Ref. 41). The 
loss coefficient is defined as: 


Ер - Dr? 
qı 


Ptl was mass averaged at the inlet. Pt2 was calcu- 
lated assuming that the static pressure vas atmos- 
pheric and adding the mass averaged dynamic pressure 
derived from the measured outlet velocity profile with 
constant density. Unfortunately, the calculations did 
not give good results since the numerator was a small 
difference between large numbers, each of which was 
approximate. Thus, the uncertainty was of the same 


order of magnitude as the loss coefficient itself. 
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B. COMPARISON WITH CODE PREDICTIONS 

Three codes have been used to calculate the viscous 
behavior through the cascade for the conditions of the 
experiment; namely, a non-interactive boundary layer code, a 
strongly interactive inviscid-viscous boundary layer code, 
and a Navier-Stokes code. 

1. Non-Interactive Boundary Layer Code 

The boundary layer growth on the suction surface at 
8) = 40° and 8, = 43.4° was calculated using a boundary 
layer code developed by McNally (Ref. 7). This was the code 
used by Sanger (Ref. 5) in designing the blade shape. 
Momentum and displacement thicknesses downstream of 
reattachment were input together with the edge velocity 
distribution that was measured by the LDV. 

The results are shown plotted in Figure 72. The 
lines represent the calculated thickness and the symbols 
represent the measured thickness. There was good agreement: 
between the calculated and measured data over the forward 
half of the blade, but further downstream the measured data 
departed significantly from the calculated thickness. At 
the trailing edge, the measured boundary layer was twice as 
thick as the calculated boundary layer. The disagreement 
may be the result of the particular turbulence model used in 
the code. Figure 73 examines this suggestion by comparing 
the calculated boundary velocity profile with the measured 


profile. The differences are seen to be significant and 


164 


BETA = 40 


0.08 0.12 


0.04 


THICKNESS (") 





E > = c 3 
"e 01 02 03 04 05 05 07 085 09 | 
X/C 
BETA = 43.4 


0.12 


THICKNESS ("’) 
0.08 


0.04 


0.00 





0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 i 
X/C 
Figure 72. Comparison with Suction Surface Boundary Layer 
Thicknesses Calculated Using BLAYER 





Displacement: Caliculateds: 9 LDV 
Momentum: =- ~ Calculated. + LDV 


165 


d° 


AGT U3I^ poanseeow *"pe3ge[no[?o a 
sp EE EE S9TiJOXd 
Азтоотәл 194e71 AÁaepunog 22e831ng UOIJDMS YITA мозтледиоОр  *c/ eanDrJ 





a/a oo 
9'O 90 O0 20 0 : a/a 


MO 
20 


yO 


B 
o 
n/A 
n/A 





166 


result from the assumption in McNally's code that the veloc- 
ity profile follows a power law.  Sanger and Shreeve (Ref. 
45) used McNally's Rm with an inviscid code and 
compared measured surface pressure distributions on the suc- 
tion surface with calculated results. The calculated data 
deviated from the measured data as one moved downstream 
towards the trailing edge. The inaccuracy apparent in the 
viscous calculation would help explain the disagreement. 
2. Interactive Inviscid-Viscous Boundary Layer Code 
An incompressible viscous-inviscid strongly inter- 
active (S.I.) code being developed by Cebeci (Ref. 16) ws 
used by Krainer (Ref. 17) to predict boundary layer distri- 
butions on the suction and pressure surfaces of the present 
blade. The displacement thickness distributions calculated 
for each of the three inlet angles on the suction and pres- 
sure surfaces are shown in Figures 74-75, respectively. The 
solid lines represent the calculated thickness and the 
symbols represent the measurements. Also, shown on the two 
figures are the locations of laminar separation (LS), tran- 
sition (T), reattachment (R), full turbulence (Tu), and 
turbulent separation (TS), as given by the code, The 
predictions using the code disagreed quantitatively and 
qualitatively with the measured data as follows. 
£ -. he design inlet angle (8) = 40%), the calculations 
indicated that there vas no laminar separation on the 


suction surface. Instead, the laminar boundary layer 
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Figure 75. Comparison with Displacement Thickness on the 


Pressure Surface Calculated Using S.I. Code. 
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undervent transition to turbulence through a boundary 
layer transition. Onset of transition was calculated 
to be at 37% chord. The transition extended to 55% 
chord, where the flow became turbulent. In contrast, 
the LDV data clearly indicated the presence of a 
separation bubble, and this was consistent with the 
flow visualization results reported by Sanger and 
Shreeve (Ref. 45). 

On the pressure surface, the calculations indicated 
that transition followed laminar separation at all 
three angles. At each angle, laminar Separation was 
predicted at 76 + 1% chord, reattachment at 852 15 
chord and the boundary layer to become fully turbulent 
at 99$ chord. In contrast, the measurements showed no 
laminar separation on the pressure surface at the 
indicated inlet angles. Transition extended to no 
more than 60% of chord where the velocity profiles 


were fully turbulent. 


c. The measured boundary layer displacement thick-ness 


everywhere on the suction surface was significantly 
larger than the calculated displacement thickness at 
each inlet angle. 

On the suction surface, at the highest inlet angle 
(81 = 46°), the calculations indicated that separates 


of the turbulent boundary layer would occur at 


ДИ. 


approximately 76% chord. In contrast, the measured 
data showed the flow to be fully attached at these 
test altera: 

Clearly, the S.I. code, as it vas used for the above 
calculations provided an inaccurate prediction and descrip- 
mon of transition. Therefore, to improve the calculated 
results, the onset of transition and the transition length 
were forced as inputs into the interactive calculations. 
The results for the boundary layer growth on the pressure 
surface are shovn plotted in Figure 76, The location of the 
onset of transition and the transition length are shovn in 
the figure. The lines represent the calculated distribu- 
tions and the symbols represent the measured data. The cal- 
culated momentum and displacement thicknesses are seen to 
agree quite vell vith the measured thicknesses. 

Forcing transition on the süúuction surface did 
not give equally good results. Transition through the lami- 
nar separation (separation bubble) could not be induced at 
the lower inlet angle (8) = 40°). Figure 77 shovs the suc- 
tion surface boundary layer thicknesses for g} = 46° with 
transition modified to obtain a reasonable match to the 
measured data close to the bubble. The calculated boundary 
layers are much thinner than were measured, and the 
calculations predicted the absence of separation, which was 
consistent with the LDV observations. However, the growth 


was not significantly changed from that calculated with the 
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BETA = 46 
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Figure 76, (cont.) Comparison with Pressure Surface 
Boundary Layer Thicknesses 
Calculated Using S.I. Code with 
Forced Transition. 





Displacement: Galeulated. 979 LDV. 
Momentum: ...— alculated?7: DV. 
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BETA = 46 
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Calculated Suction Surface Boundary Layer 
Thicknesses Calculated Using S.I. Code vith 
Forced Transition. 
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5 0 fied code (Figure 74). The disagreement betveen calcu- 
lation and measurement is seen to Start immediately after 
reattachment. Presently, there is no procedure available 
within the code which could force the boundary layer to be 
thicker by a change in the input parameters. A modification 
in the model used for calculating the bubble or other im- 
provements are required before the interactive code can 
describe the observed suction surface flow behavior 
accurately. 
3. Navier-Stokes Code 

A well-tested Navier-Stokes code for cascade flow 
255000 51) was applied by Ок. Shamroth and S. K. Choi of 
Scientific Research Associates, Inc. to compute the flow 
field at each of the three inlet angles. They were supplied 
with the cascade geometry, inlet angles, stagnation pressure 
and temperature, and outlet static pressure. A free stream 
turbulence level of 2.5% was assumed, following the far up- 
stream value given by Sanger and Shreeve (Ref. 45). The 
code was run with the choice of parameters shown in the 
upper part of Table 12. The code assumes transition onset 
on the suction surface at the laminar separation and at Reé* 
= 700 for the pressure surface as shown in the table. The 
location of transition calculated by the code, as a result 
of the criteria adopted, are shown in the lover part of 


Table 12. A Prandtl mixing length turbulence model was used. 
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Preliminary comparisons of the data vith the predic- 
cions of the code, run under the conditions that no informa- 
tion from the experiment vas provided to the operators, are 
shown in Figure 78 and Figure 79. There was good agreement 
between the measured and calculated pressure coefficient as 
shown in Figure 78 for all three inlet angles. A mismatch 
is noticed at the leading and trailing edges, although the 
mismatch is relatively small. The calculated boundary layer 
velocity profile on the suction surface at 90% chord agreed 
very well with the measured profile as shown in Figure 79, 

It is certain that an input of information from the 
experiment, the transition length for example, would improve 
the comparison shown in Figure 78, and this is planned. 
Davis (Ref. 52 and 53), Hah (Ref. 54), Rai (Ref. 55), and 
Delaney (Ref. 56) have also expressed interest in applying 
or modifying Navier-Stokes codes to analyze the present sent 
CD cascade. All codes are based on the Reynolds-averaged 
Navier-Stokes equations, and therefore require the inclusion 
of empirical models for laminar-to-turbulent transition, and 


for turbulent transport of momentum and energy. 
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BETA = 46 DEG. 
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C. 


SUMMARY ASSESSMENT 


The following is a preliminary assessment of the 


numerical codes as cascade flow predictors: 


Ae 


The non-interactive boundary layer code predicted 
boundary layer growth quite well on the forward half 
of the suction surface of the blade, but quite poorly 
on the rear half of the suction surface. The inter- 
active boundary layer code predicted boundary layer 
growth on the pressure surface of the blade quite well 
once onset of transition and transition length were 
introduced as input data. 

The interactive boundary layer code predicted the flow 
field on the suction surface quite poorly. 
The Navier-Stokes code predicted the flow field in the 


cascade quite successfully. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


LDV flov measurements vere carried out to determine the 
passage flov field, the suction and pressure surface bound- 
ary layers and the vake flov field of a controlled diffusion 
compressor cascade. The measurements were performed at nomi- 
nally constant Mach number (70.25), constant Reynolds number 
(7700000) and at three inlet angles; namely, design condi- 
tion 8] = 40°, and positive incidence inlet angie -mE TE 
43.4°, 46°. The following were concluded: 

a. A separation bubble existed on the suction surface, 
close to the leading edge. The flow underwent transi- 
tion In the free shear layer above the bubble. The 
reattachment point moved dovnstream as the inlet angle 
vas increased, 

b. The turbulent boundary layer on the suction surface 
vas fully attached throughout the inlet angle range. 
Separation took place on the blunt trailing edge 
curvature. 

c. The suction surface boundary layer thickness increased 
dramatically with inlet angle and reached 15$ chord at 
the highest inlet angle (8j = 46°). 

d. The pressure surface boundary layer changed very 
little with the inlet angle variations. The flow 


underwent natural transition at about mid-chord. 
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e. Wake velocity LDV measurements compared well vith pre- 
vious calibrated pneumatic probe measurements. The 
profile of the trailing edge suction and pressure sur- 
face boundary layers remained well defined in the wake 
profiles. Flov reversal vas measured close to the 
iEzurling edge (within 0.04 chords, or 1.6 trailing 
edge diameters, downstream the trailing edge), consis- 
tent with separation occurring not before the trailing 
edge curvature. 

The LDV data were compared with the predictions of three 
computational codes: a non-interactive boundary layer code, 
a research viscous-inviscid interactive code, and a near- 
production Navier-Stokes code , The folloving vere 
concluded: 

a. The boundary layer code predicted boundary layer 
growth properly at the forward half of the blade but 
was very inaccurate towards the trailing edge. mags 
suggested that transition was modelled well but that 
assumption of a power law Shape tor the prorile was 
too restrictive. 

D. The interactive code lacked accuracy in its modelling 
ö5“eransition. Hovever, vhen the measured onset of 
transition and transition length vere input into the 
code, the pressure surface boundary layer was 


predicted vell, 
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c. The interactive code did not predict transition and 
boundary layer grovth properly on the suction surface. 
Parameters could not be input or easily edited into 
the code to correct the prediction, as vas done on the 
pressure surface, since the bubble structure vas cal- 
culated by the code vithout reference to empirical 
inputs. 

d. The Navier-Stokes code, without reference to the re- 
sults of measurements, predicted the measured surface 
pressure distributions very well at all three inci- 
dence anqles. Suction surface boundary layer growth 
was predicted well although some disagreement vas 
found in the velocity profiles near to the trailing 
edge. The entry of experimental information into the 
code is thought certain to improve the degree of 
agreement. 

What has been clearly shown, is that currently available 
numerical codes for viscous cascade flow depend critically 
on the modelling of the transition process. What was seen 
in the experiment was that the boundary layer transition can 
be followed and mapped using the LDV technique, whereas the 
laminar bubble, with its free shear layer transition, is 
invisible to the LDV unless an auxiliary seeding technique 


is found to seed inside the bubble. 


186 


Therefore, the follovinq recommendations are made con- 
cerning the future of the immediate experimental and 
computational work: 

NE A 27777 for seeding the bubble should be sought, 
and some definition of the bubble region should be 
attempted. 

2. Data should be obtained at one higher angle, (^48?) at 
which separation should move forward, providing an 
even more difficult test case for codes. 

3. The S.I. code should be re-examined in the light of 
the present DE aman and the transition model in 
the code should be revised, 

271he inputs to the Navier-Stokes code should be 
modified, in light of the test data, and the three 
cases rerun in the attempt to obtain closer aqreement 
at each angle. Comparisons with other Navier-Stokes 
codes should also be pursued. The data included here- 
in are sufficient to justify computer experimentation 
with transition models. 

Finally, it is clear that the LDV technique can be used 
to provide the detailed data necessary to calibrate the 
empiricism inherent in viscous flow computational codes. 
More data of the type obtained in the present study are 
rather urgently required, preferrably using a wide variety 


of cascade blade shapes. Also, in a separate experiment, a 
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very much larger model should be used to create a separation 
bubble which, with the reattachemnt and downstream viscous 
development, can be resolved adequately. Until proper 
models are proven for the transition process, the computa- 
tional prediction of stall in the cascade cannot be expected 


to be accurate. 
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APPENDIX A 


BÖ 1610 05 DISTRIBUTION MEASUREMENTS 


Al. INTRODUCTION 

The oil droplet particle size distribution oroduced by 
the seeding apparatus was measured in a preliminary 
experiment. Since smaller particles follow the flow better 
than larger ones, the average particle diameter affects the 
accuracy of velocity measurements and the size distribution 
affects the accuracy of turbulence measurements. The objec- 
tives of the experiment were to determine the optimal seed- 
ing pressure and seeder length, and to determine the 
particle size distribution produced by the chosen seeding 


probe configuration. 


A2. DESCRIPTION OF THE EXPERIMENT 
A2.1. Particle Generator 

The particle qenerator which was used is shown 
17 Figure 9. Commercial olive oil vas the seeding fluid. 
The oil reservoir served both as a relaxing chamber for the 
droplets and as a sump for oil overflow from the atomizer. 
Oil is fed to the atomizer where the fluid is broken into 
fine droplets. A mixture of droplets and air flows out ot 
the atomizer, through the relaxation chamber to the seeding 


probe. 
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A2.2. Particle Sizer 
The particle sizer used in the experiment was a 
MALVERN 2600 series manufactured by MALVERN Instruments of 
Malvern, England. The principle of operation is that POE 


is scattered by a particle at an angle which is a function 


of the particle's diameter. Figure Al shows a schematic 
diagram of the system. A He-Ne laser beam passes through 
the test section where particles are introduced. The scat- 
tered light is collected by 31 concentric rings. The dis- 


tribution of light intensities over the concentric rings is 
interpreted as particle size distribution through software 
that was provided by MALVERN. A detailed description of the 
system (used to measure the particle size distribution in 
the exhaust of a solid propellant rocket motor) is given by 
Pruitt (Ret,.,T57). 
A2.3. Procedure and Program rfi sr non 

First, particle size distribution was measured 
with varying pressure drop across the atomizer and with 
various length of the exhaust tube. Data were collected as 
pressure was varied between 30 and 50 psi, for three 
different exhaust tube lengths, namely, 40, 80, and 120 
inches. A vacuum cleaner constantly sucked the air from the 
test section to clear olive oil fumes and to prevent 
residual droplets from previous measurements from affecting 
the next one. A three-minute interval vas maintained 


betveen measurements to make sure the test section vas 
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clear and background measurements vere frequently taken for 
the same purpose, 

Second, the seedinq probe vas connected to the 
exhaust tube coming out of the particle relaxing chamber. 
Pressure drop across the atomizer was varied between 20 and 
50 psi, and particle size distribution data were collected 


in the Same manner ase berore. 


A3. RESULTS AND DISCUSSION 

Data from the first part of the experiment are sum- 
marized in Table Al. The data values are the percentage of 
particles with diameters smaller than the indicated values 
(in microns) for different pressures (p) and exhaust tube 
lengths (Ly. It vas concluded that the best size distribu- 
tion was obtained for L= 120", p = 40 psi. Under those 
conditions the particle generator produced droplet diameters 
below 2.6 microns, with more than half below 1.2 microns in 
diameter. The seedina probe was, therefore, designed to 
have a total length less than 120". 

The data from the second part of the experiment are 
summarized in Tables A2 and АЗ and plotted in Figure A2. 
Table A2 contains the data output by the particle sizer data 
acquisition system. Table A3 shows a presentation of these 
data in the form of a histogram of size distribution, 
average droplet diameter, and standard deviation for each 


case. Figure A2 is a plot of the histogram using the data 
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TABLE Al 
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in Table A3, The average particle diameter varied betveen 
37 and 0.93 microns. The standard deviation varied be- 
tween 0.43 and 0.58. The smallest average particle diameter 
was achieved with low pressures (20-25 psi) but the distri- 
bution was not favorable since there was a large number of 
larger particles. Forty percent of the particles were 
larger than 1.2 microns. The average diameter was favorable 
at these pressures, but the standard deviation was 
unfavorable. In contrast, for pressures between 30-50 psi, 
little change in the average diameter and standard deviation 


was noted. 


A4. CONCLUSIONS 
a. The average particle diameter for an atomizer pres- 
sure drop of 30-50 psi was 0.9 - 0.95 microns. 
b. The standard deviation in the size at those condi- 
tions was 0.44 - 0.47 microns. 
c. Ninety percent of the droplets were smaller than 
ES microns.: 
d. The data agreed with the atomizer manufacturer's 
manual which specified an average particle diameter 
of 0.8 - 0.9 microns for most mineral and vegetable 


oils at 30-50 psi pressure drop. 
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APPENDIX B 


PRELIMINARY MEASURMENTS OF THE BOUNDARY LAYER 
ON A FLAT PLATE 


Bl. INTRODUCTION 

Measurements were made of the distribution of velocity 
and turbulence in the boundary layer on a flat plate. 
Results were compared with the Blasius solution for a lami- 
nar boundary layer on a flat plate and with published exper- 
imental results. The experiment served to verify the use of 
the LDV system for boundary layer measurements in the 


compressor cascade wind tunnel. 


B2. DESCRIPTION OF THE EXPERIMENT 
B2.1. Wind Tunnel and Model 

The arrangement of the test section and model 
are shown in Figure Bl. The wind tunnel's adjustable walls 
were set to 40 degrees and the test blades were removed from 
the test section. A specially designed flat plate model was 
installed at “the 7th blade station (1/3 tunnel width). A 
cross-section of the flat plate is shovn in Figure B2. The 
section had a sharp edqe and a blunt edqe as shovn in the 
figure. The model was black anodized to reduce light 
reflections from the surface which allowed flow measurements 
closer to the wall. The flat plate was installed with the 


sharp leading edge forward, with the flat test surface set 
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Figure B2. Flat Plate Model. 
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at an angle of 38 degrees, (The flat plate surface vas thus 
at a neqative anqle of attack of 2 degrees vith respect to 
the inlet walls as shown in Figure Bl.) This prevented a 
leading edge flow separation that could cause transition to 
turbulent flow as was found by McGuire (Ref. 43). 
poss LDV System 

The LDV system was as shown in Figure 8 except 
that the Bragg cells and the frequency shifters were 
removed. 

B2.3 Procedure and Program of Measurements 

Data were obtained for a range of Reynolds numbers 
by setting the wind tunnel blower to different speeds. At 
each speed, the LDV system was moved to the desired position 
downstream of the plate leading edge. The LDV traverse was 
adjusted to position the measuring volume at the surface, 
and the coordinates were recorded. LDV data were then 
acquired as the measuring volume was manually traversed to 
select displacements along a line normal to the plate 
surface. Five (5) different boundary layer distributions 
were measured, corresponding to different blower speeds and 


(x) position downstream of the leading edge; namely, 


Boundary layer 21 x = 1.0" Rex = 37500 
Boundary laver #2 x = 1.0" Rex = 55700 
Boundary layer £3 xa- 2.0" Rex = 79000 
Boundary layer #4 20-00: Rex = 152500 
Boundary layer #5 x = 9.0" Rex = 530000 
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The two components of the velocity and the turbu- 
lence level vere measured vith the LDV system, Plenum pres- 
sure vas measured using a water manometer. Data were then 
corrected and normalized to the free stream velocity using 


the plenum pressure. 


B3. RESULTS AND DT5CU?91 

Free stream turbulence was measured to be 1.4 + 0.2%. 
Results are shown plotted in Figures B3 and B4. The 
velocity component parallel to the flat plate, normalized by 
the edge velocity is shown plotted as a function of the dis- 
tance away from the wall in Figure B3. The lines renresent 
the theoretical Blasius solution for a laminar boundary 
layer and the symbols represent the measured values. The 
data are shown plotted in linear coordinates (not in semi- 
log coordinates) so that differences between the measured 
velocities and the theoretical distribution are clear. 

It can be seen that at lower Reynolds numbers (37500, 
65700, 79000) there was good agreement between the measure- 
ments and the theoretical analysis. At higher TE num= 
bers (152500, 530000) there was good agreement close to the 
wall (below 0.5 mm) but as the measuring volume was 
traversed away from the wall a departure occurred. At a 
Reynolds number of 152500 the flov vas apparently in 
transitfen. Wang et al (Ref. 27) shoved similar results 


obtained with a hot-wire anemometer. They observed that 


20/2 


¢ RE=530000 
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Figure B3. Boundary Layer on a Flat Plate. 
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Figure B4. 
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Boundary Layer on a Flat Plate. 
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with a free stream turbulence level of 2%, transition began 
at Reynolds numbers of 100000 to 200000. In particular, at 
REX = 107000, a lack of agreement with the Blasius profile 
was observed at the edge of the boundary layer. Blair (Ref. 
26) showed that at a free stream turbulence level of 2%, the 
flow transitioned at a displacement thickness Reynolds 
number of 800. The displacement thickness Reynolds number 
of Curve #4 in Figure B3 is 766. 

The measured turbulence level distributions are shown 
in Figure A4. The data by Wang et al (Ref 27) showed 
similar trends. The measured laminar boundary layer 
turbulence levels ranged from 5$ close to the wall to 1.4% 
at the edge of the boundary layer, whereas for the high 
Reynolds number (530000) the measured turbulence was 
somewhat higher (9%). Even higher levels of turbulence 
(11%) were reported by Wang et al (Ref. 27) for the laminar 


layer when measurements were made closer to the wall. 


Bae CONCLUSIONS 
a. Measurements of purely laminar boundary layers with 
the LDV system showed good agreement with the 
theoretical Blasius profiles. 
b. Transition of the flow was detected at Reynolds 
numbers which were consistent with previously 


reported results. 


Ce 


In view of the observed agreement, it was concluded 
that the seeding particles followed the flow. 
Using the Blasius analysis, particles in the flow 
were found to follow an acceleration of a = 6200 
п/зес 2. 

Turbulence measurements vere in agreement with 
Similar previously reported meaSurements obtained 


with a hot-wire anemometer. 
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ABPENDIX-C 


İ? HCULATION OF THE INLET FLOMN REFEPENCE VELOCITY USING 
PLENUM PRESSURE AND TEMPERATURE MEASUREMENTS 


Cl. INTRODUCTION 

Each measurement of velocity was divided by a reference 
velocity to give dimensionless velocity data. The reference 
velocity was the free stream inlet velocity to the test 
section at the time the LDV measurement was made, calculated 
using plenum pressure and temperature measurements. The 
following is a description of the method used for 


calculating the reference velocity. 


G2. METHOD 
The relationship between the upstream stagnation (pt}) 


and static pressure (pj), and the dimensionless velocity, 


p (V/¥2C, Ty). assuming isentropic perfect gas flow, is 
given by 
— E E - v(X) C(1) 
vhere 
1 
75:57:75 


Losses in the quide vanes are represented usinq the loss 
coefficient uq defined as 
Pto - Ptl 
5747 0: ии 
= Ҹо 
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so that Equation C(1) becomes, assuming no density change 
across the vanes, 
Ptl ^ PtQ - 0094 = Pro - 9160828) * ug С(2) 
Losses in the test section are represented similarly as 


_ Ptl - Pt2 
A A 


Therefore, assuming no density change across the test 


blades, 
28 
= E cos”el s 
Ptl = Pe2 * djup Dä toy EEN 
COST 2 
giving 
28 
pPI = bo to; (= 0 C(3) 


cos? 82 


Using Equation C(2), Equation C(3) and Equation Gaul 


—... 
Tə Fe T MOSS OMEN A v(X) C(4) 


Pto 7 di cos281 + Wg 


Since py 1S approximately atmospheric pressure, wp and 
Wg are less than 0.1, and the velocities to the test section 
are less than 100 m/sec (X «0.1), Equation C(4) can be 
Simplified to obtain 
p Pro <ip 
A A C(5) 


in which Cg is a constant for small varios rio opsi r 


Equation C(5) relates the dimensionless upstream 
velocity in the test section to the difference between 
plenum and atmospheric pressures, and the absolute level of 
atmospheric pressure, through an unknown constant Co. Co 
depended on the geometrical configuration of the tunnel, and 
was established by calibration for each new setting of the 
inlet angle. The calibration involved measuring vj using 
the LDV system while recording (pto ~Pa), Pa and Teo. The 


dimensionless velocity was calculated using 


51 E C(6) 
2:00 vas calculated using Equation C(5). 
Mithmeguescaolrshed, Lor each LDV measurement of 
velocity within the test section, the reference inlet 
velocity was calculated using Equation C(6) and Equation 
C(5), using plenum and atmospheric conditions at the time of 


the measurement. 


APPENDIX D 


LISTING OF SOFTWARE 


CASCADE FORTRAN 


X X % X X X X X % % X % X % % % %í X % % X X % X % X % %⁄ % X % % % % % X X % X X % o ə % X X % % % % X % X X X 


x x 
x x 
XX 
x x 
x x 


THIS FILE CONTAINES A PROGRAM THAT READS THE LDV DATA CORRECTS Ba 
THE VELOCITIES, CALCULATES THE TURBULENCE AND CALCULATES THE i 
STREAM FUNCTION FIELD FOR STREAMLINES PLOTS. STREAMLINES ARE j 
CALCULATED BY D(PSII) = V*DX ASSUMING PSII = O AT THE 
LEFT BOUNDARY. 


X X % X X X % X X X X X X X X % X % X % X X % X X € X % X X X X % X % % X X X % 36 % % % % X X X X X X X % % X % % % X % X % X X % % % X X X 


* 


QD OK OK OK OK OK OK OK OK OK OK OK OK OK OK Ok Ok 


ABA ues . FLOW ANGLE 
CP12 ...... COEFFICIENTS 
PAE T STATIC PRESSURE. 
PSIE ...... NORMALIZED STREAM FUNCTION. 
Dope m STREAM FUNCTION. 

SE PLENUM PRESSURE. d 
ROUE ZAA s LDV OPTICS ROLL ANGLE. 

E PLENUM TEMPERATURE. 
TURBe ates TURBULENCE LIVEL. | 
TURBU dee U COMPONENT TURBULENCE. 
THU RAB N SS V COMPONENT TURBULENCE. 

də s VELOCITY. 

Һали. VELOCITY. 

Ката E E COORDINATE. 

E COORDINATE. 
EE ee LDV OPTICS YAW ANGLE. 


REAL CC1000)0,PSIIC100,10020,VCC100)2,XC100,100)2,X1C1002, YC100) 
REAL CO(100,20),C1C100,20),C2(100,20),C3(i00,20),PSIECI OG 100 7 
INTEGER N(100) 
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L 


1 


READ C11,x) M,CP12 

WRITE (21,x*x) ' BLADE TO BLADE FLOW FIELD ' 

WRITE (21,х) ? i ae ee U 

WRITE (21,40) 

WEITE (21,30) M 

WRITE (21,40) 

MIENTE (21,*) ' YCINCH) XCINCH) U 117 -EC 3 7 5 EL HS 
İB EES(DEG) TURB(%)' 

WRITE (21,¥) ! 


MAIN LOOP 


15 
11 


DO 10 J 


= 1)М 


ESIICI,J) * 0.0 


READ 
ROLL 
PAMB 


(11,%) Y(J) N(J) Í PAMB , ROLL 
= ROLL*3.149167180. 
= PAMBx13.596 


DO 11 I = 1,N(J) 


CORRECT 


NX 
LC 
CALL 


READ (11,x) XCI,J0,U,TURBU,V, TURBV, PT, TT, YAW 
YAW = YAW*3.1416/180. 

FÜR Ü,/V ....... UC,VC(I) 

BETA = PT7(CPAMB+PT) 

CALL SOLVE (BETA,CP12,XX) 

VELMES = XXXSQRT(1115.5*(TT+460.,)) 

TURB = SQRT(((TURBU%U)%x2+(TURBVx*xV)%*%*2)/Z2.) 
TURB = TURB/VELMES 

UC = (UxCOS(ROLL)-V*SIN(ROLL))ZVELMES 

UC = UCxCOSCYAW) 

VC(I) = (UxSIN(ROLL)+V*COS(ROLL))ZVELMES 


VEL = SQRT(UCx%2+VC(I)x%<2) 

ALFA ATAN(UCZVC(1)) 

ALFA 226507” 1200212105? 

X1(I) = X(I,J) 

TEZTEZ) С0Т0 15 

URINE CC25 21) 56« 09959001 CIX4UGSVCCIO,VEL,ALFA,TURB 
GOTO 11 

WRITE (21,320 XT61) UC, VCCİ),VEL,ALFA,TURB 
CONTINUE 

NCJ) 

55-1 

ICSCCU (X1,VC,NXFC,IC,IER) 


"I H 


0021 La LLC 


BE-SXICITIO-XI1CIO 
His = Ces I 


mI = IC + II 
PSIICI+1,J) = C(III )xDxxa/q,4C(II)xDxx3/3.4C(1)xDxxZ/2, 
*XVCCCIOxD«-PSTTCT..1) 
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C 201) 
C261757 
CUCA 
COCI U 
12 CONTINUE 
10 CONTINUE 
WRITE (21,40) 
WRITE (21,x) " STREAM FUNCTION PRINTOUTS ' 
WRITE (21,x) " ' 
WRITE (21,40) 
WRITE (21,*x) " YCINCH) X CINCH) PSII ! 
WRITE (21,x) ' ' 
WRITE (21,40) 
ICOUNT = O 
polii EC AEN, 4 
DO I4 I-s 1,NCJO 
ICOUNT = ICOUNT + 1 
PSIE(I,J) = PSII(I,;J)ZPSIICN(OJ)53) ! 


C(III) 
0119 
Gör 

VC CI) 


C LE (T-GE 21796011u55 
WRITE (21,34) Y(J),X(I,J) PSII(I,J) 

C GOTO 14 | 
C 16 WRITE (21,33) X(I,J3),PSII(1I),J) 

14 CONTINUE 

] š CONTINUE 

WRITE (21,35) ICOUNT 
C 
WRITE (97,50) M 


DO 20 J = 1,M 
NJM1= N(J)-1 
WRITE (97,51) NJM1,Y(J) 
DO 21 I = 1l,NJM1 
WRITE (97,52) X(I,J),PSIE(I,J),C0(1I,J3),C1C(1,J) CAGA 


1 C3CI,J) 
el CONTINUE 
20 CONTINUE 


30 FORMAT (6X,13,' DATA ROWS‘ ) 
31 FORMAT (3X,F6.3,6X,F6.3,6X,F6.3,6X ,F6.3,6X,F6.3,6X,F6.2,5X,F5. 2)| 
32 FORMAT (13X),F6.3,4X>F6.3,4X>F6.3,4X>F6.3,4X>F6.2,5X»F5.2) 
33 FORMAT (13X,F6.3,4X,F6.4) 
34 FORMAT (3X,F6.3,4X>F6.3,4X»F6.6) 
35 FORMAT (€2X,'# OF DATA POINTS =',14) 
40 FORMAT (7/7) 
50 FORMAT (2Х,12) 
51 FORMAT (€2X,13,2X,F6.3) 
52 FORMAT (2(2X,F7.4),4(2X.E11.4)) 
STOP 
END 
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BLAYER1 FORTRAN 


HEHE HH HH HIE HE HEH HE HE IE HE IE IE IE HE HH HE I EH HH HI IE X NX XX XX XX X X XNXX XXX XX XX XNXX XXX X XX X X X ED 
xx THIS PROGRAM READS THE BOUNDARY LAYER LDV DATA NORMALIZES IT 
XX WITH THE FREE STREAM VELOCITY AND TRANSFORMS THE DATA TO COORD. 
xx PARALEL AND VERTICAL TO THE SURFACE. IT ALSO CALCULATES THE 

xx DISPLACEMENT AND MOMENTUM THICKNESS AND THE SHAPE FACTOR. 

HHH HHH HI KK HH HH HK KH HK KK IK KK KH HE KKK KK KKH IK HK HEH NX X NX X N X XX NX XXX X 


* 


ALFA .... FLOW ANGLE. 

7 .... CDEFFICIENT. 

i  ... DISTANCE FRUM BLADE SURFACE, 
DC ...... DISTANCE FROM BLADE SURFACE NORMALIZED BY CHORD. 
DENIA... DISPLACEMENT THICKNESS. 

H ....... SHAPE FACTOR. 

PAMB .... STATIC PRESSURE. 

7... PLENUM PRESSURE. 

uM... LDV OPTICS ROLL ANGLE. 

mü  ... MOMENTUM THICKNESS. 
ieee... PLENUM TEMPERATURE. 

iÜ . .. TUFKBULENCE LEVEL. 


IDE... U COMPONENT TURBULENCE. 

TURBV ... V COMPONENT TURBULENCE. 

Meee... VELOCITY 

En... VELOCITY PARALLEL TO BLADE SURFACE. 
UPARN ... UPAR NORMALIZED BY REFERENCE VELOCITY. 
IEEE... VELOCITY 

VVER .... VELOCITY NORMAL TO BLADE SURFACE. 
VVERN ... VVER NORMALIZED BY REFERENCE VELOCITY. 


7— ... COORDINATE. 
Y Ф Ф Ф 6 6 6 6 COORDINATE, 
YAW ..... LDV OPTICS YAW ANGLE. 


O XK X X X X X X X XK X X X X X K K K K K K K X * 


REAL D(100),)DC(100), TURB(100),UPAR(100),UPARN(100),VVER(100), 
1 VVERN(100) 


C 
CHARACTERx8 SURFCE 
2750  3.5x3.1416/180. 
C 
C READ LDV DATA 
Ç 


READ (31,x) YB,KN,PAMB,CP12.ROLL,ALFA 
READ (31,x) N1 


PAMB = PAMBX13.596 
ROLL = ROLEx3.14916/180. 
ALFA = ALFA*3.1416/180. 
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C 
C 
C 


OOO 


OOO 


OOO 


SURFCE= 


DO 10 


I 


"SUCTION " i | 
IF (K.LT.0) SURFEE = "PRESSURE® 


1,N 


READ (31,%) X»Y»U,TURBU,V,TURBV,PT>TT 


CORRECT FOR U,V 


BETA 


VELMES 
TURE CID 


UN 
VN 


UZVELHES 
V7VELMES 


UNHOR 
UNVER 
UPARCI) 
VVERCI) 


IF (K.LT.02 7/Ek EE 2750 
SQRT(Xxx248Yxx2) 


D( I 


) 


DC(I) 


10 CONTINUE 


NORMALIZE WITH LOCAL FREE STREAM DATA. 


FACTOR 
DO 11 


UPARN(I) 
VVERN(I) 


il CONTINUE 


I 


UPARCN1) 


РТИСРАМВЕРТ) 
CALL SOLVE (BETA,CP12,XX) 
= XX*SQRT(1115.5%(TT+&660.)) 


= SQRTCCCTURBUXU )*¥*2+(CTURBV¥V)*¥*2)/2. )/VELMES 


(UNYCOS(ROLL)-VNxSİNCROLL))xCOSCYAM) 
UNXSINCROLLO + VN¥TSOCROLL) 

= UNHORXSINCALFA)+UNVERXCOSCALFA) 

= UNHOR*COSCALFA)-UNVERXSINCALFA) 


Вит иб 


1,N 
UPARCIO/FACTOR 
ҸҸЕК(1)/РАСТОК 





DISPLACEMENT & MOMENTUM THICKNESS, SHAPE FACTOR | 


DELTA = 0. 
THETA = 0. 
DO 12 I = 1,N1-1 
DELTA = DELTA + (1.-CUPARNCI)+UPARNCI4+1))72. )*(DCI+1)-DCI ae 
THETA = THETA + ((1.-UPARNCI))XUPARNCI)+C(1.—-UPARNCI+1))XUPARI 
1 +1) D4 CD C1413 blr 2. 
12 CONTINUE 


H = DELTA/THETA 


PRINTOUTS 


WRITE 
WRITE 
WRITE 
WRITE 


(525707 BOUNDARY LAYER OUTPUT' i 
(5:27 * =әс=-====ед===с===с=-==-==-"*"? | 
(32,40) | 
(32,30) YB,SURFCE | 


214 x 


C () C) () 


tJ 


13 


30 
51 
Se 
> 
34 
40 


WRITE (32,40) 
WRITE (32,x) ' D/C | UPAR | VVER |UPARN |VVERN |TURB * 
WRITE (32,x) "(--) İ (--) İ (--) 1 (--) İ (--) ( (57) 
WRITE (32,x) ' | | | | | , 
DO 13 I = 1,N 
WRITE (32,31) DC(1),UPARCI),VVERCI),UPARNCI),VVERNCI),TURBCI) 
CONTINUE 
WRITE (32,40) 
WRITE (32,32) DELTA 
WRITE (32,33) THETA 
WRITE (32,34) H 
FORMAT C(15X,'Ys',F6.3,2X,A8,2X, 'SURFACE') 
FORMAT (F6.4,1X,4(F6.3,1X),F5.2) 
FORMAT (10X,'DISPLACEMENT THICKNESS=!',2X,F6.4) 
FORMAT (10X,"MOMENTUM THICKNESS=',2X,F6.4) 
FORMAT (10X,'SHAPE FACTOR=!,2X,F6.3) 
FORMAT (///) 
STOP 
END 


SUBROUTINE ANGLE CALCULATES THE BLADE SURFACE ANGLE WITH RESPECT TO 
THE VERTICAL. 


SUBROUTINE ANGLE CYB,K,ALFA) 


REAL CCG) 

кии) = 0.0 
С(2) = 0.0 
7525) = 0.0 
C(4) = 0.0 


IF (K.GT.0) THEN 
READ (1,Х) МР 
BLSE 
READ (2,x) NP 
END İF 
DO 14 I = 1,NP-1 
51 £ ST.03 THEN 


READ (1,/x) YP,C3760Z£.€1 
EESE 
PEZ2U (z:x) YE,U36- 61 
END TF 
IF (YB.GT.YP) GOTO 15 
D = YB - C(4) 
2507 1.-(3.xC(3)xDexx2:2.xC(2)xDıCC1)) 
ALFA = ATANCDYDX) 
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15 


14 


CCl 61 
CE C2 
[03053065 
CCG) = CG 
CONTINUE 
RETURN 
END 
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WAKE1 FORTRAN 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXX XX XXX XXX XXX XXX X XXX XXXXX XXX XXX XX XXX XXX Xx 
xx THIS PROGRAM READS THE WAKE LDV DATA AND NORMALIZES IT WITH 

x% THE INLET VELOCITY WHICH IS CALCULATED USING PLENUM PRESSURE AND 
Ex TEMPERATURE. FOR THE CASE OF BEAMS ROTATED BY SOME ROT ANGLE. 
ххххххххХХх ХХХ ХХХ ХХ ХХ ХХХ ХХ ХХ ХХ ХХХ ХХХ 6263 2 2636 36 3636 3636 3636 Е3636363Е3Е3Е26 363636363696963630302 


лә. 
ALFACI,J) .... 
5— — ......... 
ELEM e e e nn > 
5— EEN 
217” ......... 
7”  ........ 
2£— 7:2 ...... 


27-:0) ...... 
TURB(I,J) .... 
TURBU(I;J) one 
TURBV(I,J) ... 
57:10) ....... 
UNORMC I,J) 
7:71) ....... 
Ee Ee AJ ..... 
VELNCI,;,J) 
VELMESCI;J) 
VNORMCI,;J) 
5: ....... 
— ........... 
21710) ...... . 
ҸАМ(1,Ј)..... . 


Q) C) KX X X K XK ok OK OK OK OK OK OK OK OK OK OK OK OK OK OK K k k * O O O 


NOMENCLATURE 


HORIZONTAL COMPONENT OF THE MEASURED VELOCITY. 
ANGLE OF VELOCITY VECTOR 

VERTICAL COMPONENT OF THE MEASURED VELOCITY. 
COEFFICIENT 

# OF ROWS 

# OF DATA POINTS PER ROW 

AMBIENT PRESSURE ("H20) 

PLENUM PRESSURE ("H20) 

ANGLE OF ROTATION OF THE BEAMS. 

MEASURED STAGNATION TEMPERATURE. 

CALCULATED TURBULENCE. 

TURBULENCE MEASURED IN THE HORIZONTAL COMPONENT. 
TURBULENCE MEASURED IN THE VERTICAL COMPONENT. 
MEASURED HORIZONTAL COMPONENT OF VELOCITY 
UC(I,J) NORMALIZED BY THE INLET VELOCITY. 
MEASURED VERTICAL COMPONENT OF VELOCITY 
MEASURED VELOCITY 

VEL(T:41) NORMAETZED By INLET VELOCITY. 

FREE STREAM VELOCITY 

V(I,J) NORMALIZED BY THE INLET VELOCITY. 
HORIZONTAL COORDINATE. 

VZVT 

VERTICAL COORDINATE. 

YAW ANGLE OF LDV SYSTEM. 


REAL ALFA(150,10),PT(150,103),TT(150,10),TURB(150,10), 
1 TURBU(150,10),TURBV(150,10),U(150,10),UNORM(150,10),V(150,10) 
2 MELCTSOS 10), VEELNCISO? 10), VELMESC150,10),VNORMC150,10), 


3 Х(150,10), 


INTEGER N(10) 


OOO 


Y(10),YAW(150,10) 


READ AMBIENT PRESSURE & PRESSURE COEFFICIENT 
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READ 


(12,*) PAMB,CP12,ROT 


NCJ), YCJ) 


XCI,J),UCI,;JO),TURBUCI; JO, VCI, JO); TURBVCIS 22» 
PTCI,J),TTCI,J)I,YAWCI,J) 


ҸАМ(1,Ј)хҝ5.14616/180, 


ROT = ROT*¥3.19416/180. 
PAMB = PAMB*¥13.596 
C 
C READ LDV DATA 
E 
KEAB (12:51 H 
DO 10 J = 1,M 
READ (12,х) 
DO 11 I = 1,N(J) 
READ (12,x) 
1 
ҸАМ(1,Ј) = 
11 CONTINUE 
10 CONTINUE 
C 


C NORMALIZE AND LDV DATA 
C 


DO 127071 
DO 13 I = 1,N(J) 
BETA = PT(I,J)Z(PAMB+PT(I,J)) 
CALL SOLVE (BETA,CP12)XX) 
VELMESCI,J) - XXYSURESLMH" 5 0070160 2 
A = UCI,J)*COSCROT) - VCI,J)*SINCROT) 
B = UCI,J)*SINCROT) + VCT,J)*COSCROT ) 
A * AXCOSCYAWCI,;,J)0) 
UNORMCI;J) AVMVEENESCISJS 
VNORMCI;J) B/VELMESC I.) 
VELCI,J) z SQRT(AxxZ3BxxZ) 
VELN(I,)J) = VEL(I,J)ZVELMES(I J) 
ALFACI,J) = ATANCA/B)*180.73.1416 
TURBCI,J) * SQRTCCTURBUCI;JOX*xUCI;,JOOX*xx2*(TURBVCI,JOXVCIOMM 
1 se VEPEMESCIS JO 
TURBCI;JO = TURB(I,J)ZSQRT(2.) 
13 CONTINUE 
1 CONTINUE 


C 
C PRINTOUTS OF LDV DATA 
C 
WRITE 
WRITE 
WRITE 
WRITE 


(22,x) 
(22,x) 
(cc. 
(22,40) 


WRITE 
IB 


(677 Y 


WAKE OUTPUT' 


(LASER BEAMS ROTATED BY 45 DEG)! 


X VEEMES U-HORM V-NORM VEEN ALFA TUR 
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15 
14 
30 
31 
Sc 
33 
40 


WRITE 
1 1 

WRITE 
1 1 


(2207777119 (1N) Ch? SEC). == =< E (DEG) (5) 


Cac; )* 


mor 14 J = L,M 
21:16 —- :0) YC0)?,XC1,2) ,VELMESCİ1,4) ,UNORMC1)2),VNORMC1,2), 


VELN(1,J),ALFAC1,J),TURB(1,J) 


DO 15 I = 2,N(J) 


MWEBESCS2,319 XCISJOPVEEMESCI?JO;UNORMCI;,JO,VNORMCI?;JO), 
VELNCI,J),ALFACT,J),TURBCI,J) 
CONTINUE 


CONTINUE 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
STOP 
END 


DE eelerer LT TO TIE A 
(7X,1X,F6.3,1X,F6.l,3(1X,F6.3),1X,F6.l,1X,F6.2) 
(/6C2X,F6.1)) 

(Е9.1,4Х,РЕ6.2) 

(ZZ) 
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APPENDIX E 


ESTIMATION OF THE UNCERTAINTY IN THE MEASUREMENTS 
OF BOUNDARY LAYER THICKNESS AND SHAPE FACTOR 


Displacement thickness is defined incompressibly as 


u 


D ... (D E dr EN 


Momentum thickness is defined incompressibly as 


$ u u 
ZFZ5— —— E(2) 


It vas assumed that the uncertainty vas due mainly to an 
uncertainty of 0.002" in positioning the measuring volume on 


the blade surface. Therefore, we write 


END 
- fv u $ u > 
Sk = ei Lor nm ) dy "yal (1 - mm ии ии. 
апа 
Е(4) 
= yı 2. u Š = u u = 
s= ol (l- T 


where subscript 1 denotes the data point closest to the 
surface, 

The second integrals in Equations E(3) and E(4) are not 
changed in value by a shift in the y axis. Therefore, they 
Cannot contribute to the uncertainty in either the displace- 
ment or the momemtum thicknesses. In order to calculate the 
first term of each thickness, a linear interpolation was 


assumed for the integration. Thus, 
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u 
51 (1 - = ) ay = d 


and the uncertainty 


AĞ ox = (Lo — Ay] = 
Similarly, 

9 = fla - M ) ur dy = (1 = 

1 0 uz U 2 


and the uncertainty 


1 1 ul 
BED L 92 Su,’ 


The shape factor is defined as 


The maximum uncertainty in the shape factor is 


therefore given by 


AS ^9 
44 rm 
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